Nonabelian dark matter: models and constraints 
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Numerous experimental anomalies hint at the existence of a dark matter (DM) multiplet Xi 
with small mass splittings. We survey the simplest such models which arise from DM in the low 
representations of a new SU(2) gauge symmetry, whose gauge bosons have a small mass fj, ,$ 1 
GeV. We identify preferred parameters M x = 1 TeV, fi ~ 100 MeV, a g ~ 0.04 and the \X ~ * 4e 
annihilation channel, for explaining PAMELA, Fermi, and INTEGRAL/SPI lepton excesses, while 
remaining consistent with constraints from relic density, diffuse gamma rays and the CMB. This 
consistency is strengthened if DM annihilations occur mainly in subhalos, while excitations (relevant 
to the excited DM proposal to explain the 511 keV excess) occur in the galactic center (GC), due 
to higher velocity dispersions in the GC, induced by baryons. We derive new constraints and 
predictions which are generic to these models. Notably, decays of excited DM states x' ~ y X7 arise 
at one loop and could provide a new signal for INTEGRAL/SPI; big bang nucleosynthesis (BBN) 
constraints on the density of dark SU(2) gauge bosons imply a lower bound on the mixing parameter 
e between the SU(2) gauge bosons and photon. These considerations rule out the possibility of the 
gauge bosons that decay into e + e~ being long-lived. We study in detail models of doublet, triplet 
and quintuplet DM, showing that both normal and inverted mass hierarchies can occur, with mass 
splittings that can be parametrically smaller (e.g., O(100) keV) than the generic MeV scale of 
splittings. A systematic treatment of Z2 symmetry which insures the stability of the intermediate 
DM state is given for cases with inverted mass hierarchy, of interest for boosting the 511 keV signal 
from the excited dark matter mechanism. 

PACS numbers: 98.80.Cq, 98.70.Rc, 95.35.+d, 12.60Cn 



I. INTRODUCTION 



In the last year, it was intriguingly suggested that a va- 
riety of observed astrophysical anomalies might be tied 
together by a single theoretical framework, in which tran- 
sitions between states in a dark matter (DM) multiplet, 
mediated by new GeV-scale gauge bosons, could lead to 
production of lepton pairs [jj. These could explain ex- 
cess electron/positrons seen by the PAMELA [2], ATIC 
0, PPB-BETS [4], HEAT [5] and INTEGRAL/SPI @ 
experiments (the latter via the excited DM proposal 
(XDM) [7|]). In addition, it has been proposed that such 
transitions could account for the DAM A/LIBRA annual 
modulation via the inelastic DM mechanism (iDM) 
0. Synchrotron radiation from the leptons could ex- 
plain the WMAP haze 10]. More recently Fermi/LAT 
[TTj ] and HESS [lj] have made higher precision measure- 
ments of the e + e~ spectrum at TeV energies, confirming 
an excess above the known background, although less 
pronounced than the ATIC data. The DM explanation 
for this excess has by now been studied by numerous au- 
thors El El El EE E3, El, El, SI Si S3, SIM Si, M , 

and a plethora of models has been proposed [27j , includ- 
ing ones where the DM decays rather than annihilates 
[28(. Pulsars provide a more conventional astrophysical 
explanation 1 for many of these anomalies, but the data 



do not yet clearly prefer them over the DM hypothesis 
[3ll ]. However, constraints from secondary gamma rays 
produced by the charged leptons (or from primary neu- 
trinos) are rapidly clo sing up the allowed DM parameter 

space [pj3|M si SS si SI si, so, sn m si si 

Ea . |46|. 1471748 1 . Anticipated new data from the Fermi 
telescope is expected to tighten these constraints in the 
near future. 

The theoretical paradigm we focus on here assumes 
that the DM transforms nontrivially under a nonabelian 
gauge symmetry which is spontaneously broken below the 
10 GeV scale. Radiative corrections from virtual gauge 
bosons induce mass splittings between the DM states of 
order a g /j., where a g = g 2 /47r is the fine structure con- 
stant of the new gauge symmetry and /i ~ gv is a char- 
acteristic gauge boson mass after spontaneous symme- 
try breaking. Multiple exchanges of the light gauge (or 
Higgs) bosons gives a Sommerfeld enhancement [l|, |49| 
which can explain the large annihilation cross section 
needed in the galaxy, compared to the smaller one in 
the early universe at the DM freeze-out temperature , ex - 
pected from the relic density. In our previous paper [501 ] , 
we presented an SU(2) model along these lines which 
was designed to more easily give a large enough 511 keV 
signal as observed by INTEGRAL while also accommo- 
dating the PAMELA/ ATIC observations. 

Our goal in the present paper is to give a more compre- 



1 An even more conservative interpretation is that no new source 
is needed to fit the data; see for example ref. 12911 . or concerning 



the 511 keV excess, ref. [30 
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hensive survey of models based on SU(2) gauge symme- 
try, considering a few different possibilities for the means 
of coupling the DM to the standard model, for the rep- 
resentation of the DM multiplet, and that of the scalars 
which break the gauge symmetry. We also derive some 
new constraints on the gauge and Higgs couplings which 
are particular to this class of models. We start by dis- 
cussing a number of general issues which transcend the 
individual models. 

The paper is organized as follows. Section|TT]details the 
mechanism of kinetic mixing of dark and standard model 
(SM) gauge bosons, including its possible UV origin, and 
we derive a new constraint on the gauge coupling from 
the induced DM transition magnetic moment in the non- 
abelian case. Section HIIl discusses the alternative of com- 
munication between the dark and SM sectors by Higgs 
mixing. We derive new constraints on diagonal Yukawa 
couplings of the dark Higgs to DM, from direct detection 
and from antiproton production in the galaxy. In section 
IIVI we discuss the concept of an inverted DM mass hi- 
erarchy for boosting the predicted 511 keV INTEGRAL 
signal, and the Z 2 symmetry and nonthermal DM his- 
tory needed to make this idea work. Section [V] analyzes 
which regions of parameter space best fit the experimen- 
tal anomalies (we do not insist on explaining DAMA, 
since the constraints on the iDM mechanism have be- 
come so severe [Il|, [52j],[53j]) an d constraints from diffuse 
gamma rays, relic density, big bang nucleosynthesis, and 
laboratory constraints. 

In the remainder of the paper we discuss several spe- 
cific kinds of models, organized according to the SU(2) 
representation of the DM. Sections IVI1 IVIII and IVIIII re- 
spectively deal with DM in the doublet, triplet and quin- 
tuplet representations. In all of these models the gauge 
group is simply SU(2). For completeness and contrast, in 
section llXl we consider one model with dark gauge group 
SU(2)xU(l) and triplet DM, which illustrates the differ- 
ences between the purely nonabelian models and ones 
where gauge kinetic mixing occurs between U(l) field 
strengths. We summarize our findings in [X] Appen- 
dices [A] and [B] respectively give details of the transition 
magnetic moment and radiative mass computations, [C] 
computes the annihilation cross sections for freeze-out of 
DM in a general representation, and |E] treats the diago- 
nalization of the gauge boson and DM mass matrices for 
the SU(2)xU(l) model. 



II. KINETIC MIXING OF GAUGE BOSONS 

A simple way of generating couplings between one of 
the SU(2) gauge bosons and electrons is through non- 
renormalizable couplings of the form 



1 



4- A 



Y B^A a 



(1) 



(2) 



where A a and h are respectively triplet and doublet Higgs 
fields which are assumed to get a VEV. By having several 
triplet or doublet fields (labeled by index i) which get 
VEV's in different directions, it is possible to get mixing 
with several colors of the B gauge boson. In fl}, note that 
only a single linear combination of B vectors mixes with 
the SM. With a generic Higgs potential, we can always 
choose the linear combination of triplets in JT]) to be Ai. 
However, depending on the Higgs potential, the vector 
that mixes with the SM may be a linear combination of 
several B mass eigenstates. 

To understand the consequences of gauge boson mix- 
ing, it is useful to start with a simple example in which 
a massive abelian boson B mixes with the photon. The 
kinetic term is 

- i (F^F^ + B^B"" - 2eB^) + \i?B a B a (3) 

Since the U(l) gauge symmetry of the photon is unbro- 
ken, it must remain strictly massless. This restricts the 
form of the transformation which diagonalizes the kinetic 
term to 



eB„ 



(4) 



Therefore all particles which couple to the photon acquire 
a coupling of strength ee to the massive B gauge boson. 

For the models we consider, the mixing takes the form 
hcY^vBi" , where for concreteness we take color 1 of the 
nonabelian gauge boson to mix with the standard model 
weak hypercharge, it is straightforward to show that eq. 
P| generalizes to the similar form 



A" = A" + e cos 9 W B^ 



(5) 



where Ow is the Weinberg angle. One must further trans- 
form B\ and the Z gauge boson as 



Bl 



B± — e sin Ow 1 



Z " + ^w[-^—r 2 )B» 

m z — fi 



;i + o(£ 2 ))(6) 

;i + 0(6 2 ))(7) 



where the tilded fields are those which diagonalize the 
kinetic term. Therefore the B\ gauge boson acquires a 
coupling to the current of the Z boson, in addition to 
that of the photon. Figure [lja) shows an example of a 
X — > x'ff transition mediated by the B\. 

The mass of the Z gets shifted by a fractional amount 



Sm z 2 ■ 2 a M : 
= e sin Ow 



m z 



2m| 



(8) 



relative to its usual value. For the small values of e ~ 
10~ 3 — 10~ 4 and [i <, GeV which are of interest, this is 
a negligible shift. 
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(a) (b) 

FIG. 1: Feynman diagrams for (a) \ ~* x'ff y i a S au S e boson 
mixing (left) or (b) Higgs mixing (right). 



pressed by powers of a heavy scale, whereas models with 
SU(2) xU(l) gauge symmetry in the dark sector allow for 
renormalizable mixing of SM and dark hypercharge, in 
which case there is no reason to expect particularly small 
values of e. A phcnomcnological advantage of small e is 
that values on the order of 10~ 16 give the gauge boson 
B\ a lifetime of order 10 12 s. Such a long lifetime lets £>'s 
produced from DM annihilation propagate away from the 
galactic center before decaying. This delocalizes gamma 
rays produced by the leptonic decay products, allowing 
such models to evade HESS constraints [2(| H^|. How- 
ever we will show in section IV El that gauge bosons with 
a lifetime greater than ~ 1 s are ruled out by big bang 
nucleosynthesis for the models considered in this work. 



With gauge boson mixing, the annihilation \X ~ > 
B\B\ results in subsequent decays of B\ — ► with 
roughly equal branching ratios for all leptons I with mass 
below fi. Due to the nondiagonal couplings of Bi to the x 
states, assuming they are Majorana, there is no s-channel 
annihilation through a single virtual B\. Hence the an- 
nihilation into 4 leptons is guaranteed. For Dirac DM, 
such as in the doublet representation, this need not be 
the case, as we will discuss in section fVTl 



A. Microscopic origin of gauge kinetic mixing 

The dimension-5 operator (|TJ) can be induced at one 
loop by a heavy particle X which carries both dark SU (2) 
charge and weak hypercharge y x , if it also has a Yukawa 
coupling to the dark sector Higgs triplet. Suppose X is a 
Dirac fermion which transforms as doublet of the SU(2), 
so the Yukawa interaction is 

h x X i (T a ) i j X'A a (9) 

The diagram is shown in figure |2^a) . It generates the 
effective interaction which can be estimated as 



h x y x g 
16tt 2 M x 



Y R^A 



(10) 



so that the mixing parameter is given by e = 
h x y x gA/(16Tr 2 M x ), where A is the VEV of the triplet 
Higgs. For couplings of order unity and A ~ 10 GeV, 
Mx can be of order TeV to generate e ~ 10 -4 . 

Similarly, the dimension-6 operator @ can arise 
from a heavy doublet scalar field Si with a coupling 
A (S^TaS) (h^Tah) to another dark higgs doublet h (or 
perhaps the same one, h — > S). If S has weak hyper- 
charge y s , the analogous diagram with X replaced by S 
gives rise to the operator with A 2 = 167r 2 Mg/(<?Ays)- 



B. Long-lived dark gauge bosons 



C. Direct decay of excited DM to photon 

Because there is no mixing of By to A^ 1 in eq. ((6|) , there 
is no tree level amplitude for the decay of excited DM 
directly to a photon. For example in the case of triplet 
DM, one would have the decay %3 — * X27 if such a mixing 
existed. Instead the dominant decay is X3 ~ * X2l + l~ 
mediated by the B\. However, in the class of models 
with kinetic mixing between SM hypercharge and one 
of the dark SU(2) gauge bosons, it is inevitable for the 
single photon final state to arise at the loop level, as we 
now show. Naively, one could draw the diagram where 
form a loop connecting B\ to the photon, but this 
just renormalizes the kinetic mixing term, so it is not 
relevant. There is another process which occurs due to 
the nonabelian nature of the B±, illustrated in fig. I2fb). 

The novel feature of the gauge mixing operator is 
that B? u contains the term g{B%B% - B%Bg). There 
is thus a trilinear vertex coupling these gauge bosons to 
the weak hypercharge field strength, with strength eg. 
One consequence of this interaction is the generation of 
a transition magnetic moment for the DM. An example 
is shown in fig.^b) for the case of DM in the triplet rep- 
resentation. A magnetic moment interaction of the form 
^23X20 'tivX^F^ arises, where /i 2 3 is expected to be of 
order eg 3 / (Wir 2 M x ) . A careful computation of the loop 





It is noteworthy that pure SU(2) models generi- 
cally predict small gauge mixing parameters e, sup- 



FIG. 2: Loop diagrams which generate (a) gauge kinetic mix- 
ing (left) and (b) DM transition magnetic moment (right). 
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diagram given in appendix |A1 gives 



A<23 



eg 3 c w f ln M* 



128tt 2 M v 



(11) 



where /i is the scale of the nonabelian gauge boson masses 
and cw = cos 9w- It is straightforward to compute the 
rate I\ for xi ~* X37, 



M23 , 



r 7 = ^(<$M 23 r 

07T 



(12) 



where 5M23 = M X2 — M X3 is the energy available for the 



decay. On the other hand, the rate T 2e for xi ~* X3 e e 
is approximately (6M 23 -2m e ) 2 ((5A/23+2TO e )/(2567r 3 Mj; 
times the spin- averaged squared matrix element, 



M 



(\M\ 2 )~32gWc 2 w . 



-^(E + E_+p + -p_-m 2 e ) (13) 



(where E 1 ^ and are the energy and 3-momenta of the 
electron and positron, respectively). This varies approx- 
imately linearly over the allowed phase space, so we esti- 
mate the integral as being 

T 2e S! Ae 2 aa g {5M 23 - 2m e f(5M 23 + 2m e ) 2 /^ 4 (14) 

The branching ratio for the single photon versus the two 
lepton decay is thus 



BR-, 



c 2 w a 2 g /a H 4 (SM 23 ) 3 



8192^ 2 M 2 {5M 23 _f{5M 23+ ) 2 e/x 



ln 2 ^ (15) 



where 5M 23 ± = 6M 23 ±2m e and e = 2.71828 .... Taking 
SM 23+ = 25M 23 = 4m e but allowing for the possibility 
that SM23- -C SM 23 , we can write 



BR 7 = 2.6 x 1(T 4 



-4 a g ( A* \ 4 {1 TeVy 



100 keV\ 
' JM23- / 



200 MeV/ V M x / 



(16) 



The reference values chosen here are compatible with 
constraints which we will discuss in later sections, and 
small values of 5M 23 _ enhance the size of BR 7 . 

Even though the branching ratio for X3 ~~ > Xil due to 
the magnetic moment is small, the observable signal due 
to this process, in the diffuse gamma ray background, 
is distinctive. If the dark matter was at rest, it would 
produce a monoenergetic photon with E = SM ~ MeV. 
Since the central galactic DM has a velocity distribution 
with dispersion v/c ~ 10~ 3 , the spectrum of the photon 
is Doppler broadened with a width of order (v/c)5M ~ 1 
kcV for 8M X ~ MeV. This is just below the 1.5 keV 
resolution of SPI. The nonobservation of such a signal by 
INTEGRAL thus provides a new constraint on models 
with ^-parameter type mixing of the nonabelian gauge 
boson with weak hypercharge. 



To determine the constraint, we can compare the new 
direct photon signal with that of the 511 keV line al- 
ready observed by INTEGRAL. The latter is seen with a 
confidence level (c.l.) of 50cr and a signal to background 
ratio (S/B) of a few percent. One can predict the c.l. of 
the new signal from that of the 511 keV line through the 
relation 



(Cl.) r 



(c.l.) 5 nBR. 



(S/B) nc 
{S/B) 5U \a 



CT511 



1/2 



(17) 



where 175x1 = 5 keV is the width of the 511 keV line 
and <7 ncw = 1.5 keV is the resolution of the detector 
(which is approximately the same as the intrinsic line 
width). To understand the dependence on width, notice 
that for fixed flux, increasing the width of a line reduces 
the signal proportionally (l/c), but for fixed signal-to- 
background, it increases the counting statistics by \J~5 
since a wide line of a given intensity has more flux than 
a narrow one. These effects combine to give the \j\fij 
dependence. The background for the 511 keV line is dom- 
inated by the positronium — > 37 continuum and annihi- 
lations of positrons in the INTEGRAL telescope, effects 
which are both absent for the new signal. On the other 
hand, there is a broad instrumental line near 1.8 MeV 
which is the dominant background for the narrow galac- 
tic 26 Al line (54|, whose signal to background ratio is 
around 70/30. Putting these numbers together, and as- 
suming that (cd.)new < 3 to avoid a detection, we find 
the limit 



tg < 0.08 , 



^200 MeV \ z ( M x \ ( SM 23 _ \ 3/2 



/' 



) I 1 TeV ) I 100 keV ) 



(18) 



(recall that SM 23 _ = SM 23 — 2m e ). It is interesting that 
such reasonable values of the dark gauge coupling could 
lead to an additional signal potentially detectable by IN- 
TEGRAL. However, it would require a nonthermal DM 
history, since we will show that smaller values of a g are 
needed for the correct relic density, eq. (|36p. or in the 
case of doublet dark matter, the bound (|18p does not 
apply because the magnetic moment is suppressed by an 
additional factor of SM 23 /M X ~ 10 -6 , as we will show in 
section I VI Al 



III. MIXING THROUGH THE HIGGS SECTOR 

A. General features 

An alternative way in which the dark matter might 
couple to the standard model is through renormalizable 
operators of the form 



Xhs\H\ 2 \S\ 2 



(19) 



where H is the standard model Higgs doublet and S is a 
Higgs field which is charged under the dark SU(2) gauge 
group. If S gets a VEV vgj\p2 and also has a Yukawa 
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coupling to the DM, schematically of the form h s SxXi 
then transitions such as x ~ > xff can be mediated by 
the Higgs bosons as shown in figure [T^b). The Higgs 
sector has a mass matrix of the form 



X HS v H vs rni 



(20) 



where v H is the VEV of the SM Higgs h = y/2H. If the 
mixing is small, then the Lagrangian fields are related to 
the mass eigenstates by 



1 -t 



1 



H' 
5" 



X hs vhvs 



(21) 



Therefore S' couples with strength y/# to any SM model 
fermion / whose Yukawa coupling to H is yf. In addition, 
the H' couples to XX with strength —ysO. Thus the 
diagram involving H' exchange is of the same order in 
couplings as that with the S 1 , but at low momentum 
transfer it is suppressed by m| / mr H . 



B. Constraints on diagonal couplings 

1. No antiproton production 

An interesting qualitative difference between Higgs 
and gauge boson mixing is that in the former case, the 
Yukawa couplings are generally not off-diagonal. For ex- 
ample, triplet dark matter coupling to a quintuplet scalar 
as x a SabX h has diagonal couplings; similarly for doublet 
dark matter coupling to a triplet scalar via Xi T tjXjS a - In 
either case, the annihilation XX ~~ * S — * f + ' f~ shown in 
fig. [3Ja) occurs, resulting in quark or lepton pairs favor- 
ing the most strongly coupled fermions — the top quark. 
To avoid production of hadrons, since no antiproton ex- 
cess is observed by PAMELA, one needs to have mixing 
with a scalar that has dominantly off-diagonal couplings 
so that XiXi annihilates primarily to a pair of <S* bosons 
by virtual xi exchange. The S bosons decay nearly on 
shell and hadron production can be suppressed if the S 
is lighter than ~ 1 GeV. Note that it is impossible to 
keep the couplings strictly off-diagonal in the mass basis, 
once the relevant component of S gets a VEV, since this 
contributes an off-diagonal mass term to the DM. There- 
fore the Higgs mixing scenario in its simplest form could 
be disfavored by the lack of any antiproton excess in the 
PAMELA data. 

Moreover, diagonal couplings are constrained by direct 
dark matter searches, by the process shown in fig. [3£b). 
Translating the limit quoted in eq. (11) of ref. [l[ to the 
present case (and assuming ms = 200 MeV), a diagonal 
Yukawa coupling h s is bounded by 



0h s y N < 16tt x 10~V 



6h s <4x 10" 



(22) 



Here yiy — 10 3 is the Higgs-nucleon coupling (55|. As- 
suming that the SM Higgs mass mn is much heavier than 




N 



N 



(b) 



FIG. 3: (a) Left: DM annihilation into t-i by virtual scalar 
that mixes with Higgs. (b) Right: DM scattering on nucleon 
by scalar exchange. 



ms, this implies 



130 GeV 



2 / 1 GeV 



h s X HS < 2.7xl0~ 4 f ms V ( 
V200 MeV/ V 

To illustrate how severe (or mild) this constraint might 
be, consider the case of triplet DM Xa coupled to a 
quintuplet (traceless symmetric tensor) scalar S a b, via 
hsXaSabXb, and the cross-coupling XH 2 trS 2 to the SM 
Higgs H. Suppose for example that Si 2 gets a VEV 
S ~ 1 — 10 GeV to induce mixing with H, with mix- 
ing angle 9 = X hs Svh /m 2 H . In addition to radiatively 
generated mass splittings of the DM (as we will discuss 
below), there is a tree level contribution h s SxiX2 so 
that the mass eigenstates become linear combinations, 
X± = y/l/2(xi ± Xz)- The fluctuations SS of S12 thus 
couple to the mass eigenstates as h s 5S(x+ — X-)- There- 
fore the ground state X- can annihilate directly into a 
single fermion pair through a single intermediate scalar. 
The latter is always far off shell, so this annihilation 
channel is dominated by production of top quarks which 
hadronize and produce antiprotons, contrary to the ob- 
servations. However notice that the two diagrams in fig. 
[3ja) interfere destructively. We can estimate the effect 
of these diagrams by integrating out the intermediate 
scalar, and using the fact that mg <C itih, to get the 
effective dimension-6 operator 



0yth s 



M 4 



xxtt 



(24) 



On the other hand, the annihilation x.X ~^ S"S" by x ex_ 
change can be estimated from the dimension-5 operator 



h 2 

—xxS 12 



(25) 



Assuming that the initial x's are nonrelativistic, the ratio 
of the corresponding cross sections is of order 



o{xX -> it) & 2 y? rn% 
a( X x -> S'S') ~ h 2 M 4 



(26) 



The top quarks decay to b quarks before hadronization, 
and each b quark produces ~ 4.5 antiprotons (using 
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MicrOMEGAs [56|), so the number of antiprotons per 
positron is of the same order. The observed flux of an- 
tiprotons to electrons is approximately 10 -3 , and given 
that no antiprotons in excess of standard expectations 
are observed, we should demand that the ratio (|26|) not 
exceed this limit. For definiteness, if M = 1 TeV we 
obtain the rather weak constraint 



Afi-s 85 h s 



M 



1 TeV 



vs 



(27) 



Both (f2"T| and the direct detection constraint (f2"3")l can be 
satisfied using reasonable values of the couplings. 

Furthermore, if there are additional contributions to 
the DM mass splittings, it is possible to parametri- 
cally suppress the diagonal couplings. For example, 
consider a second quintuplet Higgs T a b with coupling 
h T Xa.T a bXb, and a VEV which splits the \ masses di- 
agonally, h T T(xiXi ~ X2X2)- In this case, the % mass 
eigenstates are not maximal mixtures of the flavor states; 
rather x+ = Xi + $X2, X- = X2~Sxi, with S = h s S/h T T 
(assuming 8 is small). If the |T| 2 |i/| 2 coupling is negli- 
gible, then the overall effect is to reduce the diagonal 
couplings by the factor S, while leaving the off-diagonal 
couplings unsuppressed. 

The constraint due to the assumed lack of production 
of antiprotons would be weakened even further if the re- 
cent claim of ref. is verified. This work questions the 
assumption that the observed antiproton background is 
actually understood in terms of physics other than dark 
matter annihilation. 



X-. 



y 5M93 -lOOkeV 
2 5M 23 ~lMeV 

8M 12 ~lMeV x 2 _ 

h X l - 

(a) 



5Mi2~100keV 



(b) 



FIG. 4: (a) Left: inverted mass hierarchy of triplet DM; (b) 
Right: normal hierarchy. 



a mixing angle is required to get the Higgs lifetime to be 
10 12 s. In section|V]it will be aruged that Higgs masses in 
the range mass ms ^ 100 MeV are the most promising 
for fitting PAMEL A/Fermi observations, such that only 
the e + e~ final state is available. Using the decay rate 
T = 2 y1m s /167T, we find that 



6 x 10 



-12 



100 MeV\ 
ms J 



1/2 



(28) 



is the required value. We will show in section fV E 31 that 
such small values are strongly excluded by constraints 
on the density of dark gauge bosons, which must decay 
before BBN. 



IV. INVERTED MASS HIERARCHY AND Z 2 
SYMMETRY 



2. No two-lepton final states 

In section IVD I we will discuss the fact that recent con- 
straints on DM annihilation from the diffuse gamma ray 
background are more severe for models in which XX 
than 41 final states (where I is a charged lepton) due to 
the harder spectrum in two-body decays. This is not an 
issue when the intermediate particle is an SU(2) gauge 
boson, since its couplings are automatically off-diagonal 
and thus two bosons must be emitted in the annihilation, 
but it might be an issue for intermediate Higgs bosons 
with diagonal couplings. However, the result (|26p can be 
directly adapted to the case of decays to a lepton pair 
instead of a top pair by substituting the lepton Yukawa 
coupling for that of the top. Even for the heaviest lepton, 
r, the result is suppressed by (y T /yt) 2 — 10~ 4 . These an- 
nihilations are thus much more rare than those with the 
it final states, and do not provide a stronger constraint 
than the one derived above, even if we only demand that 
the ratio be <C 1 rather than < 10~ 3 . 

C. Long-lived dark Higgs boson 

In order to realize the long-lived intermediate state 
proposal of ref. [2(| , it is interesting to know how small 



In the following models, a recurring theme will be 
whether it is possible to have a stable excited DM state 
which is slightly lighter than the highest excited state 
(the one that decays into leptons plus ground state) . This 
"inverted hierarchy" is shown in figure (Ha), in contrast 
to the "normal hierarchy," fig. 0|b). We proposed the 
inverted hierarchy in ref. [501 ] as a means of boosting the 
galactic 511 keV signal from excited dark matter, since 
the transition X2X2 — * X3X3 requires less energy than 
X1X1 ~~ * X3X3 and therefore benefits from a larger pro- 
portion of the DM velocity distribution. 



A. Radiative mass corrections 

Let us first review the mechanism of radiative mass 
splitting of a DM multiplet by virtual massive gauge 
bosons, through diagrams like that shown in fig. [5fa). 
Although the correction to the mass is logarithmically 
divergent, mass differences between members of the mul- 
tiplet are finite. By choosing a suitable counterterm, the 
finite part which contributes to the mass splitting can be 
defined as 

m^—ag^TtTi (29) 
3 



7 



where a g — g 2 /4n and the sum runs over all the gauge 
bosons, with mass which contribute in the inter- 
mediate state. The approximation ([29| is valid when 
fj,j <C M x . Details of the derivation are given in appendix 

m 

If Higgs mixing rather than gauge boson kinetic mixing 
is the dominant portal between the dark and SM sectors, 
it is likely that the dominant souce of mass splittings is 
the tree level contributions from the Higgs VEVs. It is 
possible however that the analogous radiative corrections 
with the intermediate Higgses, fig. [5jb) , have an impor- 
tant effect. In appendix [B] it is shown that the analogous 
formula to (1251) in this case is 



3 



(30) 



where y — y/4Tra y is the relevant Yukawa coupling for the 
Higgs multiplet in the loop, and uij is the mass individual 
components of that multiplet. 



B. Z2 symmetry 

The idea of exciting the intermediate state \2 depends 
on it being significantly populated and stable on cosmo- 
logical time scales. One possibility is for it to be abso- 
lutely stable, which should be guaranteed by some sym- 
metry. Another, which has been explored in ref. [581 ] . is 
that the state is only metastable. In section IV Fl we will 
discuss that this scenario is strongly constrained by direct 
detection considerations. In this paper we will highlight 
models that admit a discrete Z2 parity, which not only 
ensures the stability of the intermediate state, but also 
forbids transitions between it and the neighboring states, 
that could be coupled to currents of SM particles. 2 The 
absence of these transitions makes the models safe from 
the direct detection constraints. (For other references 
discussing symmetries which stabilize DM, see [H|.) 



B, 




37 



\ y 



%2 



Xi 
(b) 



FIG. 5: Examples of radiative correction leading to mass split- 
tings within DM multiplet from exchange of virtual gauge 
bosons (left) and Higgs bosons (right). 



2 Such transitions, if they exist, can always mediate decays X2 
Xi + 37, as in fig. fTJa). 
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FIG. 6: Transitions between quintuplet states mediated by 
the three gauge bosons Si, B2, B3. 



The simplest example is triplet DM \i m which only 
one gauge boson color, say B2, mixes with the SM hyper- 
charge. In this case we can assign conserved Z2 charges 
to the fields 



Xi, X3,Bi,B 3 



(31) 



and to no others. Suppose that xi is the ground state 
and X3 the heaviest state. Because of the Z2 symmetry, 
X2 can never decay into xi phis SM particles. It could 
in principle decay into X1B3, but this is kinematically 
blocked by the mass of the B3 . From the point of view of 
the symmetry, there is no light particle that can appear 
in the final state to compensate the Z 2 charge of the xi ■ 
Alternatively, we can state the condition that would 
make it impossible to keep the intermediate state stable. 
From the above argument we see that a necessary re- 
quirement is to be able to assign Z2 charge to the ground 
state. Therefore the highest excited state of interest must 
also be charged. If any gauge boson which mediates 
transitions between the intermediate state and either of 
the charged states mixes with SM hypercharge, then Z2 
charges cannot be consistently assigned. 



C. Z2 x Z 2 symmetry for quintuplet DM 

The issue of having a stable intermediate state does not 
arise for DM in the doublet representation, but it can be 
applied to higher representations, such as the symmetric 
tensor (quintuplet). We can label the canonically nor- 
malized states of Xab by 

/A-B/V3 C D \ 

(X) = C 2B/V3 F (32) 

V D F -A-B/VZJ 

The transitions mediated between these states by the 
three Bi gauge bosons are shown in figure [6l Let us 
consider how to assign Z2 charges to the states in a sys- 
tematic way. First, suppose that one of the gauge bosons, 
say B a , mixes with SM hypercharge. Then B a must not 
carry Z2 charge, while the other B's do; call these Bi. 
This implies that some subset X of x's which appear 
only linearly and not bilinearly in the gauge interactions 
of the BiS should also be charged. The states in X must 
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also have the property that they only appearly bilinearly 
and not linearly in the interactions of B a . 

Using this logic, we can make an exhaustive list of 
the possible Z2 charge assignments for a given choice of 
the B that mixes with hypercharge, which we denote by 
B a «-> e + e~. In the process, we discover that actually 
the global symmetry is larger than just Zi\ for a given 
subset X of x states, its complement Y could also have 
been chosen. This means that we can assign Z2 charge to 
states in X, and a separate Z' 2 to states in Y. Meanwhile, 
the two gauge bosons other than B a transform under 
both Z2 and Z' 2 . The result is 

B l ^e+e--. X = {C,D}, Y = {A,B,F} 
B 2 ^e+e-: X = {C,F}, Y = {A,B,D} 
B 3 <-> e+e- : X = {D, F}, Y = {A, B, C} (33) 

It turns out that the A and B states always mix to form 
the heaviest {A') and lightest (B 1 ) mass eigenstates. We 
therefore take the heaviest state relevant for the INTE- 
GRAL transition to be one in Y, and this dictates that 
the intermediate state whose stability is to be guaran- 
teed is the lightest one in X. The Z2 x Z 2 symmetry 
then insures that the ^-charged intermediate state can- 
not decay into the Zg-charged lowest state, since both 
symmetries would be violated. We will give explicit ex- 
amples in section [VIII B 11 

In order for this to work, at least one of the Z^s must 
be left unbroken by the VEVs of the Higgs fields. If there 
is only one triplet Higgs which gets a VEV to accomplish 
kinetic mixing, this presents no difficulty since then the 
Higgs components can transform in just the same way 
as the corresponding gauge fields, preserving both Za's. 
Moreover components of a quintuplet Higgs can be given 
the same charges as the corresponding DM components, 
so one Z2 can be preserved as long as VEVs appear only 
in the X or Y subsets, but not both. 

VEVs of additional doublets break all of the discrete 
symmetries, but multiple triplet VEVs can be consistent 
with the symmetries if they are orthogonal. Consider two 
triplets with VEVs A and A' in the 1 and 2 directions, 
respectively, and suppose that A is used to generate ki- 
netic mixing between B\ and the SM. Then a single Z2 
is preserved, under which the fields B2, B3, A2, A3 and 
A[ change sign, while Ai, A' 2 and A 3 do not. Adding 
a third triplet A" with VEV in the 3 direction is also 
consistent with the Z2, if A" transforms under it. 



D. Nonthermal history 

Even though Z2 symmetry guarantees the stability 
of the intermediate state \2, it cannot prevent deple- 
tion of its density in the early universe, through ex- 
actly the same process needed for the INTEGRAL signal, 
namely X2X2 -> X3X3 followed by X3 -> e+e~xi decay. 
Even more simply, the depletion could occur directly by 



X2X2 - * XiXi- I n re f- H3j we noted that this deple- 
tion could be prevented if the \' s were produced out of 
thermal equilibrium rather than through the standard 
freeze-out. If the x' s are decay products of a supermas- 
sive scalar S, their initially high energies suppress the 
annihilation cross section sufficiently long to keep the 
X2X2 ->■ X3X3 excitation or the X2X2 X1X1 relaxation 
out of thermal equilibrium in the early universe. 

In more detail, suppose that the gauge coupling a g is 
too large to yield the right relic density from freeze-out. 
It was envisioned that S could decay at a low tempera- 
ture ~ 5 MeV, resulting in mildly relativistic DM with 
momenta p ~ 10 5 T. The Sommerfcld enhancement is 
initially absent for DM with such large velocity, and in 
fact the rate of annihilations remains always less than 
the Hubble rate before cosmological structure begins to 
form, because n(av) and H both scale like T 2 . Only 
when DM begins to concentrate in halos does the rate of 
annihilations become significant. 



V. FITTING PAMELA/FERMI/HESS VERSUS 
INTEGRAL / SPI, COSMOLOGY AND 
LABORATORY BOUNDS 

A. Fits to PAMELA/Fermi/HESS 

Ref. [25| has identified regions in the parameter space 
of M x and (J ann v Ic i for the process XX — * BB (fol- 
lowed by B — + e + e~) which are compatible with the 
PAMELA/Fermi/HESS e+e - observations, as well as the 
HESS constraints on inverse Compton gamma rays pro- 
duced by the electrons and positrons coming from DM 
annihilation. 3 As we will discuss in further detail below, 
additional constraints from extragalactic diffuse gamma 
ray production favor the models in which the fl's from xx 
annihilation decay only to e + e~ and no heavier leptons. 
This implies that the mass of the B's, /x, must be less than 
twice the mass of the muon, fi < 200 MeV. The allowed 
region for this scenario is reproduced in fig. [7J It should 
be emphasized that the two-body decay XX e + e~ me- 
diated by a single B exchange is excluded because its 
electron spectrum ends too abruptly due to its near mo- 
noenergeticity [25j ; this channel also provides a very poor 
fit to the PAMELA data [11]. Moreover in the class of 
models considered here, it would be impossible to forbid 
the channels XX ~ * // where / is any SM fermion, if 
XX e + e~ is unsuppressed. 

The best fit is in the vicinity of {cr ann v Te i) = 10 -23 
cm 3 /s and M x = 1 TeV. This cross section exceeds that 
needed for the correct thermal relic density (10~ 36 • c cm 2 
[60| ) by a factor of B = 330, which is thus the required 



A recent analysis of preliminary Fermi observations of gamma 
rays from the inner galaxy is also consistent with this annihilation 
channel [47tl . 
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DM DM -* 4e, isothermal profile 
700 2000 GeV 
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FIG. 7: Best fit of ref. [13] to the XX —* 4e annihilation 
channel, in plane of {av Te i) and M x . Shaded regions in upper 
part are excluded by diffuse gamma ray constraints. 



needed to get the right relic density for triplet DM. 

Still, to avoid the stronger inverse-Compton con- 
straints on the preferred Einasto profile, it may be nec- 
essary to reduce the annihilation rate near r = 0, in 
addition to providing alternative subhalo regions for the 
annihilation. Recent work on halo formation including 
the effects of baryons indicates that the velocity profile 
steepens considerably (diverging like r -1 / 4 ) for r <, 20 
kpc instead of leveling off to smaller values [64[ as in 
pure DM simulations. (This reference also finds that the 
DM density profiles are softened near the center, a re- 
sult not corroborated by other simulations which include 
baryons [65| . but the latter work does qualitatively con- 
firm the steepening of the velocity profile [6^].) More- 
over the overall magnitude of the velocity is somewhat 
increased for r <, 100 kpc. Because the Sommerfeld en- 
hancement of the annihilation cross section scales like 
1/v, this should have a similar effect to erasing the cusp of 
the density profile, making it more similar to the isother- 
mal profile. 



B. Relic density 



boost factor, assuming a thermal origin for the DM. Even 
if a nonthermal origin is assumed, the thermal component 
should be suppressed by having an even larger cross sec- 
tion, and thus 330 should be regarded as an upper bound 
on the required boost factor. 

The example shown assumes an isothermal radial den- 
sity profile, which eases the constraints from HESS on 
the inverse Compton photons by lowering the DM den- 
sity near the galactic ridge. For preferred profiles such 
as Einasto, the fit to PAMELA/Fermi is nearly ruled 
out. The isothermal profile is considered to be unrealis- 
tically flat near the center compared to the results of the 
best N-body simulations, but it was noted in (2f| that 
long-lived intermediate bosons (the B gauge bosons in 
our case) could justify such an effective profile, due to 
the -B's traveling away from the galactic center before 
decaying [2(| • We will show that in section IV E II below 
that big bang nucleosynthesis constraints rule out such 
a long-lived B in the present class of models, hence the 
mechanism of long-lived intermediate states cannot work 
here. 

Another way of decentralizing the region of DM anni- 
hilation has been proposed in [l|, however, which could 
have a quantitatively similar effect to the softer halo pro- 
file; namely DM subhalos which populate the halo could 
dominate as annihilation sites, due to their lower veloc- 
ity dispersion and hence larger Sommerfeld enhancement. 
The small-velocity subhalo scenario has recently been 
studied in detail in ref. [45[ (see also [6l|, [H, [63[), with 
reference to models favored by the pre-Fermi analysis of 
[H[, in particular with M x = 1 TeV, fx = 200 MeV and 
a g = 0.04. This happens to be close to the preferred 
values mentioned above; we will show in the next section 
that this value of a g is just slightly larger than the one 



We have computed the early-universe annihilation 
cross section of DM in any SU(2) representation into dark 
gauge bosons. For the three representations we focus on 
in this paper, the result is 

TO 2 f 0.14, doublet ) 
Knniw) = -rrf x { 0.88, triplet \ (34) 
M x [ 5.18, quintuplet J 

Details are given in appendix [Cj Using the standard 
value (av Ye i)_= 10~ 36 cm 2 ■ c needed for thermal relic 
abundance [601 ] . comparison with the cross section |34|) 
indicates that the values of a g required are 

( 0.077, doublet ] / m \ 
a g = { 0.031, triplet } x — (35) 
[ 0.013, quintuplet J V 1 lev / 

(relic density value) 

As mentioned above, there are motivations to question 
the assumption that DM has a thermal origin, such as 
our inverted mass hierarchy proposal [5(| (see also [58T ] 
and (67jV It is important to notice that to justify a non- 
thermal origin, the thermal contribution must be smaller 
than usual so that it is subdominant to the nonthermal 
contribution; thus the annihilation cross section would 
be larger. The values (f3"fj)) should then be regarded as 
lower bounds. 

These bounds can be evaded if the DM has stronger 
Yukawa couplings to dark Higgs fields; for example triplet 
DM can have the coupling hxi^ijXj to a quintuplet Higgs 
E. If h 3> g then the freeze-out density is determined by 
h and the gauge coupling can be smaller than in (I36p . In 
such a case, it should be kept in mind that the annihi- 
lation in the galaxy will probably also be dominated by 
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Higgs boson exchange; notice that the mass scale of the 
Higgs bosons cannot naturally exceed that of the gauge 
bosons by a large factor, since the scale of spontaneous 
breaking of the dark SU(2) gauge symmetry is dictated 
by the mass scales in the Higgs sector. Thus late-time 
annihilations would likely be dominated by Sommerfeld- 
enhanced Higgs exchange diagrams. The expected boost 
factor would thus still be ~ 300 even in cases where a g is 
much smaller than indicated in (|36|) . 



C. Mass splittings and the XDM (iDM) mechanism 

In contrast to the above values of a g , the paradigm 
of ref. [l| would at first seem to suggest smaller values 
ot g ~ 10~ 3 , because the radiative mass splittings of the 
DM multiplets go like a g fj, (where ^ is the scale of the 
gauge bosons masses) and it was presumed that [i ~ 1 
GeV as the largest value compatible with no produc- 
tion of antiprotons by the decays of the gauge bosons 
after XX BB annihilation in the galaxy. Since the 
XDM hypothesis requires x mass splittings of order MeV, 
a =MeV/GeV~ 10~ 3 would be indicated. 

However we have argued above that lighter gauge bo- 
son masses fj, ~ 100 MeV are in better agreement with 
gamma ray constraints. The generic estimate 5M = 
\a.gH gives ~ 2 MeV for such masses and the preferred 
coupling a g = 0.04 from section IV Al This is in just the 
right range for having excited DM states which can decay 
to e + e~ and the ground state. 

For other applications, like the iDM mechanism for 
DAMA, or our inverted mass hierarchy variant of XDM 
[50] |. it is desirable to have splittings which are perhaps 
smaller than the MeV scale. In section IVII B 21 we will 
show that with sufficiently complicated Higgs sectors 
(three triplets in this example) it is possible to reduce 
the mass splittings below the generic level of afi. It is 
also possible to design the gauge symmetry breaking (by 
appropriate choices of VEV's or the DM representation) 
so that no x mass splittings are induced by gauge bo- 
son radiative corrections; for doublet dark matter this 
is true regardless of the Higgs respresentations. In that 
case, the splittings must come from Yukawa couplings 
and then it is possible to decouple the scale of the gauge 
boson masses from that of the splitting. 

It should be emphasized that getting the excited dark 
matter (XDM) mechanism to produce a large enough sig- 
nal to explain the INTEGRAL/SPI observations is not 
as easy as just having the right DM mass splitting; one 
must generically saturate partial wave unitarity bounds 
for the excitation cross section to get a large enough rate 
[H|. We leave the details of reanalyzing this problem 
to work in progress [6^|. The same can be said (even 
more so) of the iDM mechanism for DAMA. The region 
of parameter space consistent with the DAMA annual 
modulation as well as other direct detection experiments 
is essentially excluded [13], [Hj]. We give less emphasis 
to trying to implement the iDM mechanism. 



D. Overcoming diffuse gamma ray and CMB 
constraints 

We have already seen that constraints from gamma 
rays originating as brehmsstrahlung or inverse Compton 
scattering of the emitted leptons can often rule out mod- 
els which would have provided good fits to the PAMELA 
and Fermi observations 32]- [45]. Not only annihilations 
within our own galaxy provide such constraints, but the 
accumulated effect from early redshifts and other halos 
on the CMB and diffuse gamma ray background can be 
severe. For example, ref. [42| obtains the 95% c.l. CMB 
bound 

(crw rc i) < 4 - 8 x 10" 24 cm 3 /s (36) 

for the model with XX ~ ¥ 4e and M x = 1 — 2 TeV (where 
their efficiency factor / for transfering energy to the in- 
tergalactic medium is approximately 0.9). This is barely 
compatible with the fit to PAMELA/Fermi/HESS for the 
same model in ref. f25j , reproduced in fig. [7] 

Many papers which place gamma ray constraints on 
annihilating DM assume that only two leptons are pro- 
duced, instead of the four which are predicted by the 
class of models we are considering. Given that the pre- 
ferred models are near the borderline of being excluded, 
subject to large astrophysical uncertainties, the distinc- 
tion between the relatively hard, monoenergetic input 
spectrum for two-lepton annihilations versus the softer 
four-body final states is important. In particular, ref. 
[HI (see section 4.1.3) has quantitatively shown this to 
be the case. 

Furthermore, in excluding a given model, one should 
keep in mind the correlation between the best fit model 
parameters (the DM mass, annihilation cross section, and 
gauge boson decay branching ratios) with the assumed 
DM galactic density profile, since varying the latter can 
cause significant changes in the former. For example 
some papers refer to best-fit models as determined by 
ref. [25|, but use different DM profiles to compute the 
constraints than those used to fit the PAMEL A/Fermi 
data, making it unclear which models are really ruled 
out. 

E. Relic dark gauge (or Higgs) bosons and big 
bang nucleosythesis 

In this section we consider cosmological constraints on 
the lightest stable or metastable particle in the dark sec- 
tor. Since we have identified the mass scale [i = 100 
MeV for the portal boson as being favored by fits to the 
PAMELA/Fermi/HESS data, we will take this to be the 
lightest particle, be it the gauge boson in the case of 
gauge kinetic mixing, or a Higgs boson in the case of 
Higgs mixing. By this assumption we avoid the introduc- 
tion of any scales which are even lower than 100 MeV. 

Some of the bounds we derive implicitly assume that 
the dark gauge bosons were in equilibrium with the rest 
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of the plasma at a high temperature, so that their abun- 
dance is known around the time when they are becoming 
nonrelativistic. Even if the mixing parameter e is too 
small for interactions with electrons to achieve thermal 
equilibrium with the dark sector, one should remember 
that kinetic mixing arises from some higher scale physics, 
such as a heavy X particle which transforms under both 
the dark and the SM gauge symmetries; recall eq. ©. 
Even for small values of e, such an origin for the kinetic 
mixing can insure equilibrium between the dark and SM 
sectors at the TeV scale. 



1. Long-lived gauge bosons 

In previous sections, it was noted that dark gauge 
bosons with long ~ 10 12 s lifetimes could have provided 
an escape from gamma ray constraints on annihilating 
DM through the mechanism of ref. [2(| , but we now ar- 
gue these would also dominate the energy density of the 
universe at the time of BBN, assuming the DM was pro- 
duced thermally. Let us consider the least dangerous 
case of fi = 10 MeV gauge bosons. Further, suppose that 
the SM becomes supersymmetric above the weak scale, 
so that the number of degrees of freedom is doubled; if 
instead there is a desert of no new states, this will only 
make the BBN constraint stronger. When the DM parti- 
cles freeze out between T = M x and M x /20, they transfer 
their entropy to the dark gauge bosons. This increases 
the energy density of the latter by at most a factor of 
two, since there are more gauge degrees of freedom than 
DM ones. In the meantime, between temperatures of 1 
TeV and fi = 10 MeV, the SM degrees of freedom are 
differentially heated relative to the dark gauge bosons 
by a factor of approximately (214/11/2) 1 / 3 = (9.7) 1/3 , 
due to the change in the number of degrees of freedom 
from 214 to 11, and the fact that the gauge bosons had 
been heated by a factor of ~2 by the DM annihilations. 
(The precise value depends on the dimension cLr of the 
DM representation, but for the small-d^ models we con- 
sider, this has no effect on the ensuing bound.) Thus at 
T = 10 MeV, the energy density in dark gauge bosons 
is suppressed by a factor of (9.7) 4 / 3 = 21 per degree of 
freedom. By T = 1 MeV this suppression has gone down 
to a factor of 2.1 due to the gauge bosons being nonrel- 
ativistic. However there are 3 colors and 3 polarizations, 
so this counts as approximately 4.5 extra species, and is 
ruled out. We conclude that the gauge boson lifetime 
should be less than 1 s (the time corresponding to T = 1 
MeV), requiring that 



e > 4 x 10 



_ u /100 MeV\ 



1/2 



(37) 



We used the decay rate T = ^ae 2 fi for B — > e+e~. 

Even if the thermal relic DM density is highly depleted 
by having a large annihilation cross section, the above 
arguments hold, since most of the energy of the original 
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FIG. 8: Annihilations which deplete density of stable B2 
bosons. 



thermal DM population is deposited in the gauge bosons, 
regardless of how much DM is left. The only obvious 
way to avoid the above constraint on e is to somehow di- 
lute the original DM even more relative to the SM, e.g., 
by having even more extra degrees of freedom present 
at a TeV than in the minimal supersymmetric standard 
model. We note that the B bosons will not equilibrate 
with the SM for values of e lower than (f37f , so equilibra- 
tion cannot serve to dilute the dark gauge bosons. 



2. Stable gauge bosons 

Typically only one color of the dark gauge bosons 
mixes with the SM, say B\, while transitions between 
B2 and B3 can be mediated by the nonabelian mixing 
interaction ge cos OwF^B^B^ which we referred to pre- 
viously in section III CI leaving the lighter of these two 
states stable against decay. We must verify that its relic 
density is not too large. 

For defmiteness, suppose the stable gauge boson is B 2 . 
The most efficient process for depleting B 2 is the scat- 
tering B 2 B 2 — > B\B\, shown in fig. followed by the 
decays B\ — > e + e~. We will show that this is true even 
if B\ is heavier than B 2 . 

The cross section for B 2 B 2 — ► B\B\ can be estimated 

as 



e -AE/T ^£ I ^ 

A* 2 V M 



Sfj, 



1/2 



(38) 



where AE is the energy barrier: AE = if fi 2 > Mi> 
and AE = 26 /1 = 2(fj, 2 — fit) if fit > fJL 2 . The factor of 
(Sn/fi) 1 / 2 arises from the velocity of the final state parti- 
cles, which is 1 in the more familiar case of annihilation 
to light final states. The freeze-out temperature for this 
reaction is determined as usual by setting n b 2 (<tv) equal 
to the Hubble rate, using the equilibrium density of a 
massive particle for n b 2 i one finds that 
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(39) 



for fi = 100 MeV and a g — 0.04. This implicit equation 
quickly converges to a solution by iteration. Values of Xf 
as a function of AE/ fi are shown in figure [HI 
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FIG. 9: Freeze-out value Xf versus AE/fj, (solution of eq. 
([39)) ') for the process B2B2 — > where AE = max(2( / ui — 

M2),0). 



As long as the interactions of fig. [5] are in equilibrium, 
the abundance Yb 2 tracks that of Yg 1 , whose principal 
connection with the SM is through the decays and in- 
verse decays B\ «-> e + e~. The decay rate is suppressed 
by e 2 , and for the small values of e we obtain in the ensu- 
ing bound, it is consistent to neglect scattering processes 
B\B\ «-> e + e~ whose rate goes like e 4 . In appendix [Dl 
we show that the processes BBB — > BB are able to keep 
the gauge bosons in kinetic equilibrium with themselves 
down to a temperature given by x k = n/T^ = 17.5, so 
B\ would maintain the equilibrium abundance of a non- 
relativistic particle Y cq until this temperature. At lower 
temperatures, it disappears due to its decays: 



Y, 



y eq (^)e- r ' = r cq (^) e -( r /^))- 



(40) 



where T = ^ae 2 /j, is the decay rate, x — /i/T, and H(/j.) 
is the Hubble rate at T = \x. Since Y eq (x) 

-6 



we 



find that Y eq (x k ) ~2x l(T b . The analysis of ref. [70j 
shows that a good estimate of the relic abundance of B2 
is obtained by evaluating Yb 1 (which is the source for Yb 2 
in the Boltzmann equation) at xf. Yb 2 (00) = Yb 1 {x/). 
On the other hand, the present abundance of stable Bi 
bosons must not exceed the observed DM abundance. 
Using baryons as a reference, 



„ / \ ^dm tun 
Yb 2 {oo) < — —rib 



3 x 10" 



(41) 



where £Idm /&b — 5, r?6 — 6 X 10 , tun is the mass of 
the nucleon, and we took fi = 100 MeV. Putting these 
results together, we obtain the bound 



e > — 



1 / 3(ln(^ x 10 8 )-x k + %\nx k )) M 

x f \ 
9 x 10~ 9 



M r . 



1/2 



(42) 



Since Xf < 28.7 (the value when AE = 0), this is approx- 
imately an order of magnitude stronger than the bound 



3. Long-lived Higgs bosons 

We now consider the case where the Higgs boson S 
that mixes with the SM is the lightest metastable state of 
the dark sector. The gauge bosons provide no more con- 
straint in this case since they are presumed to be heavier, 
and although they are stable, they efficiently annihilate 
into dark sector Higgses with a negligible relic density, 
~ (/i/M x ) 2 smaller than the closure density. 

If the coupling of the Higgs to the SM is too strongly 
suppressed by the small mixing angle 9, there will be 
similar problem as the one involving metastable gauge 
bosons, discussed above. The Higgs should decay before 
nucleosynthesis to avoid dominating the energy density 
of the universe. We can directly adapt the result (|5T)) by 
replacing e — > 9, 11 — + mj, e — > y e (the electron Yukawa 
coupling, y e = y t m e /m t ~3x 10~ 6 ): 



> 4 x 10" 



100 MeV 

ms 



1/2 



(43) 



Of course, this also forbids the possibility of a long-lived 
intermediate state (2 Of for transporting them outside the 
galactic center before decaying into e + e~. 



F. Long-lived intermediate DM states and direct 
detection constraint 



In section IIVI we discussed the implications of an ab- 
solutely stable intermediate DM state, protected by a 
discrete symmetry. This symmetry also made the mod- 
els safe from downward transitions X2 — > Xi mediated by 
nuclear recoil in direct detection experiments, since the 
gauge boson B3 was forbidden from mixing with the SM. 
However, if the symmetry is not present and B% does mix 
with hypercharge, interesting constraints can arise, since 
the state \i generically has a lifetime longer than the age 
of the universe, has a significant relic density, and can un- 
dergo X2 — * Xi m the detector (53[. The latter process is 
not kinematically suppressed since it is exothermic, and 
it leads to strong constraints on the mixing parameter 
e. Ref. HI finds the 90% c.l. limit e < 2 x 10~ 6 from 
CDMS for M x = 1 TeV, SM 12 = 100 keV for the small 
splitting which would be relevant for the iDM explana- 
tion of DAM A, and (j, = 1 GeV. As explained above, we 



If Xf < Xk, then the bound is slightly modified since B\ main- 
tains equilibrium density until Xf. 



4 X 10- 9 , 
e > (17.3 
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1/2 
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FIG. 10: Experimental constraints on e versus gauge boson 
mass fx taken from ref. [7lj . Enclosed regions are excluded by 
anomalous magnetic moments, beam dump experiments and 
supernovae. 



prefer fj, = 100 MeV, which makes the constraint even 
more severe, 



e < 2 x 10" 



/' 



100 MeV 



(44) 



since the %-nucleon cross section scales like e 2 /// 4 . 

Notice that the window between (|4~4"|) and our BBN 
or relic density bounds (|37I42|) is only a few orders of 
magnitude. This region of parameter space is also below 
those which could be probed by complementary experi- 
ments, as illustrated in fig. [Hjl taken from ref. [7l| (see 
also ref. [z3|,[z3|.) In the models we consider, the bound 
(|44|) can be evaded if we insist upon the Z2 symmetry 
which forbids the transitions leading to direct detection. 
This makes it possible to have models which could also be 
probed by laboratory experiments such as beam dumps. 
Another way to evade (|4~4"|) can arise if the mass splitting 
between the intermediate and ground state is too large 
[73| . since direct detection experiments do not look for 
very large recoil energies. The inverted mass hierarchy 
could thus be useful for this purpose even if there is no Z2 
symmetry and the intermediate state is only metastable. 



VI. DOUBLET DARK MATTER 

We now begin our investigation of more specific classes 
of models, organized according to the SU(2) representa- 
tion under which the DM transforms. If the DM is in the 
doublet representation, it must be vector-like (Dirac) in 
order to have a bare mass term, 



(45) 



In this case, DM number becomes conserved. Its abun- 
dance could be due to its chemical potential rather than 



freeze-out, similar to the baryon asymmetry, and so a 
nonthermal origin could be considered more natural than 
for Majorana DM. 

There is no way to split the masses of the doublet 
through radiative corrections from the gauge bosons, be- 
cause each member of the doublet has equal-strength in- 
teractions with all three gauge bosons. For example sup- 
pose only B\ were to get a mass fii; the contribution to 



the x mass matrix is 8Mn_ 



1 ^1 



ij jk 



afiiS, 



k ■ 



But we can get a splitting through the VEV of a triplet 
via the Yukawa interaction 



hx^ a x&s 



(46) 



The suffix on A3 is a mnemonic for the fact that (for 
convenience) we take its VEV to be in the a = 3 direc- 
tion, since this gives the mass splitting 1L/1A3 between 
the Dirac states xi an d X2- 



A. Gauge kinetic mixing 

Let us first consider the case of gauge kinetic mixing 
as the portal to the SM. With the above mass splitting, 
either B\ or B2 must mix with the SM hypercharge so 
that transitions between xi an d Xi can occur, with the 
production of e + e~. The triplet VEV which generates 
the mass splitting is not suitable for generating the ki- 
netic mixing of the gauge boson via jY^Bft" A3. In 
fact, such mixing is dangerous from the standpoint of 
constraints from direct DM searches, since it would in- 
duce diagonal couplings via B% of the DM to nuclei. One 
possibility is to have an additional triplet, A", coupling 
as in eq. (JTJ) , which gets a VEV along the 1 (or 2) direc- 
tion. The extra triplet VEV serves another purpose, by 
completely breaking the SU(2) gauge symmetry, whereas 
a single triplet would break SU(2)— >U(1). Assuming that 
Ai gets its VEV along the 1 direction, the spectrum of 
the gauge bosons is 



Hi=gA 3 , n 2 = g^Ai + Aj, /j l3 =gA 1 (47) 

With this spectrum and the couplings described above, 
-B3 is stable, but B2 can decay into B3 7 via the non- 
abelian gauge mixing interaction egcwF^vBi^B^. If we 
assume that Ai > A3, then annihilations B3B3 — ► B1B1 
effectively deplete any potentially dangerous B3 relic den- 
sity in the early universe. 

Alternatively, kinetic mixing could be accomplished 
by a Higgs doublet as in eq. with VEV h — 

(v/V2)(l,l) T . This would cause mixing of only B\ to 
the SM. The gauge boson mass spectrum in this case is 



Mi = ^2 = g\ A 



= gv. 



(48) 



Similarly to the case of two triplets, B2 can decay into 
B37, but for this spectrum, the B3B3 — > B\B\ annihi- 
lation channel is kinematically blocked. Therefore the 
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FIG. 11: Spectrum of doublet \ states (left) and two possi- 
bilities for gauge boson spectra discussed in the text (center 
and right). 





FIG. 12: Canceling loop diagrams contributing to the transi- 
tion magnetic moment of doublet DM. 



mixing parameter e must satisfy (|42|) to effectively de- 
plete the relic B^s. The mass levels for the DM and 
gauge boson states are summarized in figure 1111 

As mentioned in the previous section, doublet DM has 
the advantage of allowing the gauge coupling to be as 
large as needed for getting the right annihilation cross 
section, without additional constraints from the size of 
the X1-X2 mass splitting, SM = 2/1A3. For example one 
can adopt close to the preferred value from section [Vj 
a g = 0.054 (notice that this gives the correct relic density 



for doublet DM, eq. (36|i. if M x = 700 GeV, compatible 
with the allowed region in fig. [7]), and take the gauge 
boson mass at the 100 MeV scale, assuming the argument 
of the previous section that diffuse gamma ray constraints 
prefer the 4e annihilation channel over an admixture of 
e and fi. The triplet VEV's are then of order /j,/g ~ // ~ 
100 MeV, and the Yukawa coupling should be h ~ 10~ 2 
to accommodate the excited DM (XDM) mechanism for 
explaining INTEGRAL/SPI. 

A distinctive feature of the doublet DM model is that 
its transition magnetic moment is suppressed relative to 
that of triplet DM. The diagrams which contribute are 
shown in figure [TSJ The group theory factors from the 
DM gauge couplings of the two diagrams are respec- 
tively T22T21 and T2\t"x\ (the Pauli matrices), which are 
equal and opposite. Therefore the sum of the diagrams is 
suppressed by the small mass difference 5M/M and the 
gamma ray line from \i ~> Xi7 decays is too weak to be 
detected by INTEGRAL. The constraint a g < 0.08, eq. 
(TTS)) . does not apply. 



B. Higgs mixing 

Since we have argued that a DM Yukawa coupling 
to a triplet is already necessary to get the doublet mass 
splitting, it is tempting to make a more economical model 
without gauge kinetic mixing, by letting this triplet mix 
with the SM Higgs through a A|#| 2 |A 3 | 2 coupling. The 
most stringent of the constraints on A from section Ull Bl 
is (|2U|) . arising from direct detection of the DM. For h — 
10" 2 this gives A < 3 x 10~ 2 if m A3 = 200 MeV and 
(A3) = 1 GeV. Saturating this inequality leads to the 
mixing ang le 6 = 2 x 10~ 4 , according to (|22), This value 
is consistent with our BBN constraint (|4"3"|) . 



C. Diagonal couplings to Bz 

Even though we took care to avoid the direct anni- 
hilation channel X1X1 ~* e + e~ through virtual B3 pro- 
duction, by forbidding mixing between B3 and the SM 
hypercharge, it is impossible to forbid X1X1 ~ * BiBz- If 
only B\ couples to the SM but not B2, this results in the 
final state e + e~ B2, where B2 is invisible. In the foregoing 
we have noted that the two-body final state e + e~ is ruled 
out, because its spectrum has the wrong shape to fit the 
PAMELA and Fermi observations. The three-body final 
state is much more similar to the four-body one in this re- 
spect, however, because the two visible leptons share the 
energy of the incoming x's with the i?2- They thus have 
a soft spectrum which is qualitatively similar to that of 
the four-body case. Thus the XX ~ * e + e" B2 channel in 
this model is on a similar footing to the XX ~^ 4e one in 
models with Majorana DM. If B2 also mixes with the SM 
hypercharge so that B2 — > e + e~~ , they become identical. 



VII. TRIPLET DARK MATTER 

We now take Xa to be a real (Majorana) triplet of 
SU(2). It can have a bare Majorana mass Mx a X a - I n 
this case, mass splittings can be generated radiatively, as 
well as at tree level. A doublet VEV h gives equal con- 
tributions to all the gauge boson masses, so it does not 
generate any mass splittings between the Xa' s - It is thus 
more economical to assume there is at least one triplet 
VEV contributing to the SU(2) breaking. However this 
is not enough to fully split the DM states, since a single 
triplet VEV would leave two of the gauge bosons degen- 
erate in mass, and the radiative corrections would then 
do likewise for the DM states. We are led to introduce 
either a second Higgs triplet as in the doublet DM case, 
or a quintuplet. It is also interesting to consider a model 
with three triplet VEVs, since this gives additional free- 
dom in arranging the DM spectrum to have an inverted 
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or normal hierarchy. 5 In the following, we consider these 
different Higgs sectors and the \ a mass splittings that 
arise due to radiative corrections. 



A. Two triplet Higgs fields 

Let us turn on VEV's for two triplets in orthogonal 
directions, A} and A|, for example. It is easy to write a 
Higgs potential whose minima have this property: 
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V = A,(A^ - vlY + A 12 (Ai • A 2 ) 2 (49) 



As long as A12 > 0, the energy is minimized for orthogo- 
nal VEVs. 



Mass spectra 



The gauge boson mass spectrum is 



fii=gA 2 , H2=gAi, fi 3 =gJA{ 



(50) 



(With no loss of generality, one can take A;[,A 2 > 
by doing a global gauge transformation.) The radiative 
corrections to the DM masses are 



SM 2 
5M 3 



ga I Ai 
ga I A 2 



1 

2 
1 

2 

--<?a(Ai + A 2 ) 



A? 



A| 



Al 



A| 



(51) 



Depending on the ratio A2/A1, this can correspond to 
either the normal or inverted hierarchy. To see the range 
of possibilities, define Ai = Acos# and A 2 = Asin#, 
and subtract from each 5Mi the average splitting (since 
this just renormalizes the bare value M x ): 



5Mi 
SM 2 
SM 3 



gaA(l + cos9 -2sin9) 
gaA (1 + sin 9 - 2cos6>) 
#oA (-2 + cost? + sin 6») (52) 



The spectrum is plotted as a function of 9 (which lies 
in the range [0, tt/2] due to our requirement that both 
VEVs be positive) in figure [T21 There it is clear that 
the inverted hierarchy occurs if 6 <C 1 (A 2 <C Ai) or 
9 = ir/2 (Ai <C A 2 ), while the normal one occurs if 
9 = 71/4 (Ai = A 2 ). 



5 We prefer the inverted hierarchy since it can boost the effective- 
ness of XDM, and avoid the direct detection constraint (I44[ i. 



FIG. 13: Mass splittings of triplet \ states for the model 
with two orthogonal triplet VEVs as a function of 9 = 
tan _1 (A 2 /Ai), eq. (J52J) , in units of ag^J A\ + Af. 



2. Inverted hierarchy and Z2 symmetry 

To discuss the phenomenology of this model, we must 
specify which of the gauge kinetic mixing operators 
Ai ■ B^F^" I 'Ai are assumed to be turned on. The sim- 
plest possibility, and the one that allows for Z 2 symmetry, 
is that only one of them is significant, say the one corre- 
sponding to Ai. Then only B\ mixes with the SM, and 
we can assign Z 2 charges to B 2 , B 3 , X2, X3- The un- 
charged state xi cannot decay into \2, so to implement 
the inverted hierarchy for INTEGRAL, we should choose 
9 < 7r/2 to make \i the intermediate state. 

By choosing 9 = tt/2, hence Ai -C A 2 , we obtain from 
(fS"0"|) the gauge boson mass spectrum fi 2 < fj,\ < [i 3 . Ac- 
cording to the argument of section fVE 21 the gauge mix- 
ing parameter must then exceed the lower bound (|42l) . 

To be compatible with a nonthermal origin, the gauge 
coupling must be larger than a g = 0.03, according to eq. 
([36]). Taking a g = 0.06 and (j, ^ gA = 100 MeV, for ex- 
ample, the gauge coupling is then g = ^ATia g = 0.87, 
and A = fx/g = 115 MeV. The largest mass split- 
ting is of order |a/i = 3 MeV. We have the freedom 
to adjust the smaller splitting as desired by choosing 
9 = tan- 1 (A 2 /Ai). Taking 9 = 0.4tt gives SM 23 = 2.1 
MeV and 6M13 = 150 keV, which is small enough to com- 
fortably enhance the 511 keV signal to the level observed 
by INTEGRAL [11]. 



3. Normal hierarchy 

Since it might be argued that the window for iDM 
to explain the DAM A/LIBRA annual modulation is not 
completely closed [74j, for completeness we consider the 
case of the normal mass hierarchy. As is clear from fig. 
IT3l it arises from choosing 9 close to 7r/4. If one wants 
to have both the iDM and XDM effects for DAMA and 
INTEGRAL, respectively, then B 3 must mix with the 
SM hypercharge, in addition to B\ (if 9 > 7r/4) or B 2 (if 
9 < 7r/4). Interestingly, the spectrum ([4"8"|) shows that the 
boson which does not mix with the SM is the lightest one, 
while B 3 is the heaviest. Therefore the lightest gauge 
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boson is stable and the BBN bound (|42f on e applies to 
this model. 

With the normal hierarchy there is no requirement 
for a nonthermal origin of the DM, so we consider the 
value a g — 0.03, eq. (|56"|) , needed for the correct ther- 
mal relic density, and the boost factor 300 needed for 
PAMELA /Fermi. The gauge coupling is g = 0.61. Fig. 
1 of ref. [l| shows that this value of the coupling gives 
approximately the required value of the boost factor 
for jU < 1 GeV, assuming the DM velocity dispersion of 
a = 150 km/s. We are free to adjust the triplet VEVs 
to obtain the desired mass splittings. For example with 
SM 13 ~ SM 23 = 2 MeV, one finds A x - A 2 = 750 MeV. 
To get a small mass splitting 8Mx2 = j9 a g{^-i ~ A 2 ) of 
order 100 keV, if one wishes to explain DAMA, the two 
VEVs Aj and A 2 have to be tuned to be equal to each 
other to within one part in 70. 



(covering 1/8 of a sphere), one determining the orienta- 
tion of A 2 (covering 1/2 of a circle), and two controlling 
the orientation of A 3 (on a full sphere). Depending on 
the vacuum state, the DM mass splittings can take on 
any hierarchy. Unlike the previous case, \3 need not be 
the heaviest DM state. 

At a generic position in the vacuum space, all the 
B mass eigenstates mix with the SM hypercharge vec- 
tor. Assuming mi < m 2 <C TO3 in a normal hierarchy, 
one can implement both the iDM and XDM dark mat- 
ter excitation mechanisms. There are transition mag- 
netic moments between all pairs of Xi and Xj allowing 
for single-photon decays of \2 and X3- The exception 
is for Aj = A3 = 0; in that case, there is an unbroken 
Z2 symmetry, and only B\ mixes with the SM. Then ei- 
ther iDM or XDM is possible (not both), and one of the 
excited DM states will be stable. 



4- Nonorthogonal VEVs 



2. Orthogonal VEVs 



For a generic Higgs potential, Ai and A 2 are not or- 
thogonal {e.g., A 2 could be nonzero in both the 1 and 2 
directions.). Aside from changing the details of the gauge 
masses and fermion mass splittings, this has the same ef- 
fect as turning on kinetic mixing terms for both B\ and 
_B 2 . This is because the mass eigenstates of the vectors 
become mixtures of these two directions; the gauge inter- 
actions of the mass eigenstates with the \a can be put in 
canonical form with a corresponding rotation of xi,2- In 
this basis, the \ a are mass eigenstates. However, the B 
vector that mixes with the SM hypercharge vector is a lin- 
ear combination of both the B\ and i? 2 mass eigenstates, 
so both the \i <-* X3 and Xi X3 transitions couple 
to the SM. This implies there is no Z 2 symmetry pro- 
tecting the intermediate state, in the inverted hierarchy 
case. Thus it is important to keep the VEVs orthogonal 
in that case, whereas relaxing this assumption does not 
hurt the normal hierarchy scenario. In fact it makes it 
simpler, by requiring only a single gauge kinetic mixing 
term to be nonnegligible, while still coupling both DM 
transitions to the SM, as required by the iDM and XDM 
mechanisms. 



B. Three triplet Higgs fields 

1. Generic VEVs 

For a generic Higgs potential with three triplet Higgs 
fields, we can use gauge transformations to align A\ in 
the 1 direction, with A} > 0, and A 2 in the 1-2 plane 
with A 2 > 0, while the direction of A3 remains gen- 
eral. The vacuum manifold can thus be parametrized by 
an overall amplitude A = (J2 i lA^ 2 ) 1 / 2 and five angles: 
two measuring the relative amplitudes of the three VEVs 



The potential (HHJ) can be generalized to one that leads 
to three orthogonal triplet VEVs: 



V = ]>>(A 2 -, 2 ) 2 +L(A 



A3) 2 



eye. perm. 



(53) 

As long as the rji couplings are positive, the desired vac- 
uum state is a minimum of the potential. There is no 
obstacle to assuming that only B\ mixes with the SM 
vectors, if one wants to incorporate the Z 2 symmetry 
that prevents Xi from decaying. We will show that both 
normal and inverted DM mass hierarchies are possible, 
with xi as the intermediate state. Many configurations 
of the VEVs are compatible with the iDM and/or XDM 
mechanisms. 

Taking each of the respective fields Aj to align along 
the ith direction, the gauge boson masses are given by 
/ii = g\J A\ + A| and cyclic permutations, while the x 
mass splittings are 8M\ = — ka g (p2 + Ma) P ms cyclic 
permutations. The full range of possibilities for the spec- 
trum can be explored by parametrizing the VEVs in 
spherical coordinates, 

Ai = Asin(0) cos(4>) 
A 2 = Asin(0) sin(0) 
A 3 = Acos((9) (54) 

where A = (X^^?) 1 ^ 2 ' an d the angles are restricted 
by < 8 < tt/2 and < < tt/2 so that each A l is 
positive. The resulting SMi's, shifted to set the average 
value i J^i $ M{ to zero, are shown in figure [T4l 

To obtain the inverted hierarchy, where Xi is th e in- 
termediate state with a mass close to the heaviest state, 
one possibility is to take A3 -C Ai, A 2 , which essentially 
reproduces the two Higgs case studied above. This corre- 
sponds to = 7r/2 and <f> < tt/2 in fig.[IH with X2 and X3 
being respectively the ground state and highest excited 
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FIG. 14: Spectrum of triplet \ states for model with three 
triplet VEVs parametrized as in (|54|l and with the average 
mass shift subtracted. SMi are in units of ga g A/2. Each 
panel shows a different value of 6. 



state. The other possibility is to take A3 3> Ai > A2, 
corresponding to 9 = and cj) < 7r/4, in which case X3 is 
the lowest mass state. From fig. [TJ] one can also see ex- 
amples of the normal hierarchy, for example near 9 = ir/2 
and <f> — 7r/4. 

It is interesting to notice that smaller mass splittings 
than the generic scale a g /i can be obtained near spe- 
cial values of the VEVs. When Ai and A2 are equal, 
<p = it/ 4, the masses Mi and M2 become accidentally 
degnerate. By tuning the VEVs to be close to this point, 
the 100 keV scale desired for the iDM splitting can be 
achieved even if a g \i has the right magnitude for getting 
the XDM splitting. 



C. Quintuplet Higgs field 

Allowing DM to couple to a quintuplet Higgs field, 
which is a symmetric traceless tensor E Q b, gives further 
flexibility in model building, since the pattern of gauge 
boson masses induced by the E VEV is different than 
for triplets, and one also has the possibility of a Yukawa 
coupling yXa^abXb t° give tree-level contributions to the 
DM mass splitting. In addition one still wants at least 
one triplet or doublet Higgs to generate kinetic mixing of 
one of the -B's to the photon. 

For a fairly general class of potentials, the VEV of E 
can be chosen to be along the diagonal components. Let 
us take this as a simplifying assumption and show how 
much can be accomplished with just the two components, 
which we denote by 



A — B 



2B 



-A-B 



(55) 



(note the different normalization of B than for the corre- 
sponding quintuplet \ field in (|32[) ). If (E) is much larger 
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FIG. 15: Mass splittings of triplet \ states for model with 
two components of quintuplet VEV, eq. (|55[) . parametrized 



by tan 6 = B/A. 5Mi are in units of \ga g ^J A 2 + B 2 

than the triplet or doublet VEV, the resulting spectrum 
of gauge boson masses is approximately 

IH=g\A + 3B\, fX2=g\2A\, M3 = <? |A - 35| (56) 

where we used the generators 



Sac) 



(57) 



where "~" denotes that the expression must be sym- 
metrized on ab and ce. The radiative mass corrections of 
the DM states are 

SM 3 = ~ga g (\A + W\ + \2A\) 
SM 2 = -lga g (\A + 3B\ + \A-W\) 



-ga g (\A-3B\ + \2A\) 



(58) 



Parametrizing the VEVs by tan 9 = B /A, we obtain the 
full range of possibilities by letting 9 range from to 7r. 
Although the region ir/2 < 9 < it can be mapped onto 
0<#<7r/2by gauge transformations, the freedom to 
do so is generally inhibited by VEVs of other Higgs fields 
such as triplets, which we mention below. The result is 
shown in fig. 1151 



1. Normal hierarchy 

Fig. 1151 shows that the normal hierarchy occurs in the 
region of 9 = and 9 = 7r/4 (with the same shape of 
spectra at 9 = tt and 9 = 37r/4). For illustration con- 
sider the case of 9 = 0, which corresponds to B < A. 
The gauge bosons have the spectrum /13 <, /ii < /i2- The 
order of the DM masses is M\ < A/3 < M2, and we can 
turn on transitions between Xi~X2 and Xi~X3 which cou- 
ple to the electron vector current, by mixing B3 and B2 
with SM hypercharge. This can be accomplished using 
two triplets, which we will call A2 and A3, with VEVs 
A2 and A3. A model-building challenge is to find a scalar 
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potential which gives rise to this symmetry breaking pat- 
tern together with that assumed for the quintuplet. 

We noted above that it is easy to make a potential 
for triplets that gives rise to orthogonal VEV's. Sup- 
pose we do this; then global SU(2) transformations can 
be used to orient them in the 2 and 3 directions, re- 
spectively. Next consider the £ sector. The term 
A(tr£ 2 — v 2 ) 2 is 0(5) symmetric under rotations of the 
vector (A,B/y/3,C,D,F), where C,D,F are the off- 
diagonal components of S. To break this symmetry in 
such a way as to prefer the A, B components, we can add 
terms 

A 2 AfSA 2 +A 3 A^SA 3 (59) 

which are linear in A and B when the Aj get their ex- 
pected VEVs, and thus lead to nonzero VEVs for A 
and B. This would be spoiled by a term of the form 
A^EA 3 , but the latter can be forbidden by separate dis- 
crete symmetries under which A 2 or A 3 change sign. 
These symmetries are weakly broken by the gauge ki- 
netic mixing terms, which would presumably give rise to 
a small interaction through loops. This would 

generate perturbations to the previous analysis due to 
the presence of small off-diagonal VEVs in E a &. 



2. Inverted hierarchy 

Fig. [T5] also reveals the inverted hierarchy at 9 = OAtt, 
0.97T and 7r/2 ± e. The latter occurs when \A\ < \B\. 
Consider the case tt/2 — e where M 2 < Mi < M 3 . We 
need B\ to mix with the SM in this case, suggesting a 
triplet Ai with VEV in the component A}. One can use 
an analogous potential to ([59]), AiAf EAi + A 2 A^I]A 2 , 
to generate VEVs in the A, B components, if we add the 
additional triplet A 2 , which however does not play any 
role in the gauge kinetic mixing. 

For this scenario, the gauge bosons have the spectrum 
M2 < M3 < Mi- Thus B\ which mixes with the SM 
is the heaviest. The relic gauge boson constraint (T4"2"l) 
then applies, making this model susceptible to labora- 
tory searches for light gauge bosons that mix with the 
SM. 



VIII. QUINTUPLET DARK MATTER 

As the highest DM representation we will consider 
here, we turn to the quintuplet case, where \ab is a trace- 
less symmetric tensor. The gauge generators in this rep- 
resentation are given in (|5T[) , which for conciseness is not 
symmetrized in its indices, but the actual generator must 
be symmetrized in ab and ce, with accompanying factor 
of 1/4. We will label the canonically normalized states 



of Xab by 

(A-B/y/3 C D \ 

(X)=\ C 2B/V3 F _ (60) 

V D F -A-B/VSJ 

(Notice the change in normalization of B compared to 
our choice for quintuplet Higgs fields in {55}.) These are 
the mass eigenstates at tree level. 

A. Radiative mass corrections 

In previous sections we have given explicit expressions 
for the gauge boson masses assuming various patterns of 
symmetry breaking. Here we will leave them unspecified 
and study the x radiative mass splittings as a function of 
general values of /Xi. The x mass splitting term is given 
by 

8V mass = ~-a ^ X ab Tab t ce T LfgX Sg (61) 
d 

For general values of the gauge boson masses, we find 

SV mass = -afuKA + V^B) 2 + C 2 + D 2 + AF 2 ) 

—a^ 2 (AA 2 + C 2 + AD 2 + F 2 ) (62) 
-afj, 3 ((A - V3B) 2 + AC 2 + D 2 + F 2 ) 

In the simpler case where \i\ — ^t 3 , the terms which mix 
A and B cancel and the mass terms are diagonal. The 
average mass splitting in this case is — A[L\ — 2/i 2 . Sub- 
tracting away this central value, we obtain the hierarchy 
of mass splittings 

(A, D) : +28, (C, F) : -8, B : -25 (63) 

where 8 — 4a g (/ii — /i 2 ). In the more general case where 
fii = /i + 8i/Aa g , [i2 = fJ,, fJ-3 = M + ^/Aag, we get 
splittings equal to 

(A / ,B'):±2^S 2 -S 1 8 3 + S 2 

C:6 1 -25 s , D:<5i+<5 3 , F : S 3 - 25i (64) 

These can be parametrized using 

8i = 5 cos 8, S3 = S sin 9 (65) 

so that S 2 = 8a^[(/i 1 — fi 2 ) 2 + (a*3 — M2) 2 ] controls the 
overall magnitude of the splittings, whereas 9 controls 
the relative values (and can be in the range [0, 2tt] since 
there is no restriction on the signs of the gauge boson 
mass differences). We get 
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(66) 
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state, giving rise to D D — > CC followed by C — > B'e + e~ . 
Both of these examples have analogous counterparts just 
on the other side of the degeneracy between states. At 
8 >, 7r/4, the roles of C,F and B\,B% are interchanged, 
while for 8 ^ ir, the roles of C,D and B2,Bs are inter- 
changed. 

From fig.[T5]we identify several other possible inverted 
hierarchy realizations: 6 = 0, with C, D as the topmost 
relevant states, 8 = tt/2, involving D. F and 6* = 57r/8, 
involving C, F, and = 37r/2 with In short, near 

every place where two mass eigenvalues cross at an angle, 
one can have an inverted mass hierarchy. There are six 
such values of 8 where this occurs. 



FIG. 16: Spectrum of quintuplet states as a function of 9 
which parametrizes gauge boson splittings. Units of SMi are 
272a[( Ml - M2 ) 2 + ( At3 -M2) 2 ] 1/2 . 

The spectrum as a function of 8 is shown in fig. 1161 The 
complete range of possibilities for splitting by radiative 
corrections alone is thus encompassed in the figure, as- 
suming given values of the \ii can be achieved by the 
appropriate choice of scalar VEV's. 

B. Gauge interactions and mass hierarchies 

The interactions of the gauge bosons with the quin- 
tuplet states are off-diagonal, as expected for Majorana 
particles. Suppressing Lorentz indices and gamma ma- 
trices, they are proportional to 

Asaugo ~ B\{DC + V3FB + FA) + B 2 (2AD - FC) 
+B 3 (CA + V3BC + FD) (67) 

The transitions which can be mediated are shown in fig- 
ure El This diagram is useful for determining what kinds 
of DM mass spectra can be consistent with explaining 
the various experimental observations. 

1. Inverted hierarchy and Z2 symmetry 

Here we give some examples of the inverted mass hier- 
archy with Z2 symmetry which can help boost the pro- 
duction of low-energy positrons as observed by INTE- 
GRAL/SPI. In section IIV CI we identified the possible 
discrete symmetries which could exist for a given choice 
of gauge kinetic mixing. Consider the case where B$ 
mixes with hypercharge. According to the arguments in 
section HV C[ either D or F should be chosen as the in- 
termediate state. Fig [TBI shows that at 6 < tt/AF can be 
the stable intermediate state with a small mass gap be- 
low C . Scattering processes FF — > CC in the galaxy can 
be enhanced, followed by C — > B'e + e~ via B3 exchange. 
Similarly, at 8 tt, D can be chosen as the intermediate 



2. Normal hierarchy 

If one prefers a model with normal mass hierarchy, 
hg- EH shows that there are several possibilities, close 
to points where the B' curve is tangent to that of C, D 
or F. These occur at 6 = 0, 7r/2, 5tt/4. Notice that very 
small mass splittings can be arranged near these points 
with relatively little tuning of 8 due to the fact that the 
curves are tangent to each other. This gives another way 
of obtaining smaller splittings than the generic size. 

Curiously, in no case can the heaviest state A' be rel- 
evant for XDM, because there is no gauge interaction 
which couples it to the lightest state B'. Instead A' is a 
spectator, and the highest relevant state is either C, D, F, 
one of which happens to be degenerate with A' at the an- 
gles 8 = 7r/4, 7r or 37r/2. 



3. No combined hierarchy 

Because of the extra complexity of the quintuplet spec- 
trum and gauge couplings, it is tempting to look for a 
situation where both the normal and the inverted hier- 
archies could exist simultaneously, combining the advan- 
tages of the latter for XDM while still leaving open an 
iDM explanation for DAMA. This turns out to be im- 
possible, however. First consider the situation where 
only one gauge boson mixes with the SM, say B\. To 
be compatible with both iDM and XDM, B\ would have 
to mediate transitions between the ground state B', an 
admixture of A and B, and two other states. Perusal of 
the transitions in figure [5] shows that no gauge boson has 
this property. 

The next alternative is that there are two gauge bosons 
which mix with the SM, one for the iDM transition and 
one for the XDM. The problem here is that then Z2 sym- 
metry is broken for two gauge bosons. The B1B2B3 
gauge interaction then forces it to be broken for all of 
them, and no Z2 exists to protect the higher intermedi- 
ate state for XDM from decaying. 



20 



IX. AN SU(2)xU(l) MODEL 

A. Motivation 

Most of the models described above have the advan- 
tage of allowing for the inverse hierarchy of mass split- 
tings which can enhance the galactic 511 keV signal seen 
by INTEGRAL; however this comes at the expense of a 
nonthermal history for the DM in order to keep the inter- 
mediate mass state from being depopulated in the early 
universe. Furthermore, purely SU(2) DM models do not 
allow for the excitation X1X1 ~^ X2X3 which would have 
half the energy requirement of X1X1 ~~ * X3X3- I n the for- 
mer case, one need only produce a single e + e~ pair (if 
M 2 is only slightly above Mi), while in the latter, there 
must be at least enough energy for two pairs, and the 
excitation rate is therefore suppressed by the lack of suf- 
ficiently energetic DM particles in the galactic center. On 
the other hand, models with an extra U(l) in the dark 
gauge sector can have X1X1 ~~ * X2X3 by virtue of mixing 
between the gauge groups when they are spontaneously 
broken. 

Even the simplest SU(2)xU(l) model is considerably 
more complicated than most of the pure SU(2) examples. 
First, the DM is necessarily vector-like (Dirac), in order 
to have a large bare mass while carrying the extra U(l) 
charge, but the Dirac states must be split into Majorana 
states by the Higgs which spontaneously breaks the U(l). 
Furthermore, the gauge group must be completely bro- 
ken, unlike the standard model where SU(2)xU(l) breaks 
to U(l). Following [75[ we refer to this extra require- 
ment as "charge breaking." Custodial symmetry needs 
to also be broken in order for the excited DM states to 
be able to decay into SM particles, since otherwise the 
gauge bosons can be paired up into charged states such 
as = \J\j2(B\ ±i£?2), analogous to the W bosons of 
the SM. This charge is conserved if custodial symmetry 
is unbroken [the "charge breaking" mentioned above only 
insures that there is no unbroken U(l)], which would pre- 
vent the transitions between similarly charged \ states 
needed by the XDM and iDM mechanisms. Here we will 
analyze in some detail a model of triplet SU(2)xU(l) 
DM with these necessary properties, which was outlined 
in ref. 75]. The potential needed for getting the desired 
pattern of Higgs VEVs is presented there. 

B. Specification of the model 

Consider two Weyl triplets Xi an d Xi which have equal 
and opposite dark hypercharge ±y'/2. They can be given 
the bare mass term MxiXi- Once the gauge symmetry is 
broken, mass splittings can arise both through radiative 
corrections and through Yukawa couplings to Higgs fields 
which acquire VEVs. We will assume that the radiative 
corrections dominate the mass splittings which respect 
X — x' number conservation (the Dirac mass terms), while 
the Yukawa couplings (see (|69p below) are responsible for 



splitting the degenerate Dirac states into Majorana ones. 

To compute the radiative corrections, we must first find 
the spectrum of gauge bosons. As shown in ref. [ill , com- 
plete breaking of SU(2)xU(l) requires two Higgs dou- 
blets with equal and opposite dark hypercharges ±y, 
whose VEVs take the form 

/cosa\ , / 0\ 

hi = v x . , h 2 = v 2 L 68 
\ sm a J \ 1 ) 

For convenience we will also include a triplet Higgs field 
Aj with VEV (Aj) = ASn. This breaks the custodial 
symmetry at tree level. Without the triplet, custodial 
symmetry breaking first appears at one loop. To avoid 
the extra effort of computing loop effects, we parametrize 
the symmetry breaking using the triplet VEV. Finally, it 
is necessary to include a Higgs field <f> with dark hyper- 
charge —y' which can split the Dirac components of the 
DM states, so that there are no diagonal DM couplings 
of the U(l) which mixes with SM hypercharge; such cou- 
plings are strongly constrained by direct DM searches, as 
discussed in section IIII Bl The Yukawa couplings which 
accomplish this are 

hfaiXi + h'<P*XiXi (69) 

C. Mass eigenstates 

The VEVs of the four Higgs fields hi,A,cj> give rise to 
the gauge boson mass matrix in the basis B\, B2, -B3, Y 

/AO gys 2a v 2 1 \ 

A+S I 

I A + 5 gy{c 2a v 2 +v 2 ) I 

\gys2av 2 gy{c 2a v 2 + v%) B J 

(70) 

where A = g 2 v 2 , B = y 2 v 2 + y' 2 4> 2 , v 2 = vf + v 2 , 
C2a = cos 2a, S2a = sin 2a, S — g 2 A 2 and tf> represents 
the VEV of the U(l)-breaking Higgs field in (M). In 
order to give analytic expressions, we will consider the 
off-diagonal charge-breaking elements to be small pertur- 
bations ej, and the custodial breaking to be even smaller, 
6<^/\A-B\.^ 

The diagonalization of (|70[) is worked out in appendix 
IE1 The masses eigenvalues of the four gauge bosons are 
given in eqs. (|E10[) and the mixings in eq. (|E11[) . The 
standard model couplings induced by the mixing term 
eY^vYgM and the rotation matrix (jElip thus involve the 
mass eigenstates B[, B' 3 , Y', while B2 remains uncoupled 
to the SM currents. As a result, all the gauge bosons 
except B2 have relative short lifetimes due to the decay 
into e + e~. We come back to the decays of B2 below. 

The aforementioned gauge interactions give rise to 
radiative corrections to the Dirac masses of the form 
J2iXi3MiXi- Only the SU(2) gauge interaction vertices 
contribute to the splittings, because the hypercharge in- 
teractions give equal contributions to each 8 Mi. Relating 
the flavor eigenstates Bi in terms of the mass eigenstates 
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B' a by Bi = Ri a B' a , and denoting B' 4 = Y', the contribu- 
tions to the 5Mi are 



SMx = 
8M 2 = 

SM 3 = 
The coefficients R 



i 

\ a 9 ^M2 + X)^liMiJ 



(71) 

are given in eq. (|E11|) . Ignoring 
terms of 0(5 2 ), and subtracting the SM 3 contribution 
from all SMi (since we are only interested in mass differ- 
ences) we obtain 



SMi 



5Mo 



A 



A 



e 2 f(AB) 



2^ 2 VA-e 2 f(A,B] 



SMx = 



(72) 



where e = y/ej + ef, tan0 = ei/e 2 , f(A,B) = 2y/B(A- 
B)- 2 +[y/A(A - fl)]- 1 , and V = c e s e (A - B)<5/e 2 . One 
can show that the function /(A, £?) is positive for all 
values of A, B. Since we have assumed that 5 <C e 2 /| A — 
.B|, this gives a normal hierarchy with Mi ~ M2 < M3. 
Whether xi or \i is the lightest state depends on the 
ip 2 term in SM 2 . It is of order S 2 /e 4 , which can compete 
with the order S term since only S/e 2 need be small for 
the consistency of our approximations. 

The above mass splittings refer to Dirac states in the 
absence of the Majorana masses induced by (|69|) due to 
the VEV of 4>. The effect of the latter is to split the 
Dirac mass eigenvalues by ±Mh + h')(tf>). Thus we get 
two sets of states whose mass splittings are given by (1721 , 
but they are offset from each other by m = (h + h' )(</>). It 
would be consistent with our approximations to consider 
m ^S> 5 Mi, so that the more massive set of states does 
not play any role at late times. 



D. Phenomenology 

It is interesting that the mass splittings ([72")) are para- 
metrically suppressed to smaller values than the generic 
u g \x estimate, and that the hierarchy between custodial 
symmetry breaking and charge breaking, S/e 2 , translates 
into the hierarchy of masses Mi ~ M2 < M 3 , i.e., it ex- 
plains why one mass splitting should be parametrically 
smaller than the other. In terms of the model parameters, 
the smallness of et can be arranged by assuming y <C g, 
i.e., the dark hypercharge of the Higgs doublets is much 
smaller than the SU(2) gauge coupling. The even smaller 
breaking of custodial symmetry seems to require an un- 
naturally small triplet VEV A in our implementation, 
but as mentioned above, this is only a crutch to avoid 




Xi 




FIG. 17: (a) Left: diagram for B 2 
diagram for \2 — * Xi + 37 decay. 



(b) 



4e decay, (b) Right: 



loop computations, since custodial symmetry is violated 
by hypercharge interactions even without the triplet. Be- 
cause of the loop suppression on top of the small y cou- 
pling, this effect is indeed expected to be smaller than 
the charge breaking. 

Another striking point is that B 2 is the only nonabelian 
gauge boson which has no mixing with Y, due to the 
form of the mass matrix (|70[) . hence neither does Bi mix 
with the SM hypercharge (however B 2 does decay into 
electrons via the process shown in fig. [T7T a)). This means 
the transition between xi and X3 1S the only one that does 
not couple to the electron current, while the X2 Xi 
and xi X3 transitions do. The xi ^ X2 transition 
can thus be used for the iDM mechanism. Moreover, 
assuming xi is the lightest state, X1X1 ~* X3X2 (possible 
because of B3-Y mixing) followed by X3 ~^ X2^ + e~ can 
realize the XDM scenario, if SM 32 = M 3 - M 2 > 2m e . 
The other possibility is that X2 is the lowest state, so 
that X2X2 X1X3 (enabled by B3-B1 mixing) is the 
excitation channel. 

In the case where Xi is the ground state, we have 
the situation where X2 is a long-lived intermediate state 
whose transitions X2 Xi,3 couple to currents involv- 
ing nucleons or electrons. This situation is highly con- 
strained by direct detection experiments, as we explained 
in section IVFl Let us first show that X2 is cosmologically 
long-lived. The X2 state can decay to Xi + 37 through 
the operator ea 2 m~ A B 3 F^ which is induced through a 
virtual electron loop, shown in fig. [TTTb). Ref. [53[ (see 
also [7q ) estimated the rate to be 



r(x2 -> xi + 37) = 



l7e 2 a g a i {SM 32 ) 
2 7 3 6 5 3 tt 3 mf ^ 
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For the parameter values we have favored, e 
M x = 1 TeV, n = 100 MeV, SM 32 = 100 keV, a s 



(73) 

10~ 3 , 
= 0.04, 

the lifetime is around 10 21 s, much greater than the age of 
the universe. The constraint (|44|) thus applies. According 
to ref. [53] , such small couplings are inconsistent with an 
iDM explanantion of the DAMA observations. Thus the 
main potential advantage of this model over a simpler 
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one, such as doublet DM, is the lower energy threshold 
excitation channel XiXi ~> X2X3 f° r XDM. 



X. CONCLUSIONS 

We have surveyed a range of experimental and cosmo- 
logical constraints on the simplest models of dark matter 
with a hidden SU(2) gauge symmetry, with a view to- 
ward explaining the PAMELA, Fermi/LAT and INTE- 
GRAL/SPI electron/positron excesses by DM annihila- 
tion or excitation. Although new constraints on inverse 
Compton gamma rays associated with the e + e~ produc- 
tion are making it more difficult to accommodate the 
scenario, models like those we pinpoint where M x = 1 
TeV, a g ~ 0.04, the mass of the intermediate gauge or 
Higgs boson is <, 100 MeV, and the annihilations pro- 
ceed via XX ~ * 4e rather than 2e or any combination of 
heavier leptons, seem to still be viable. 

There are two uncertainties in the properties of DM ha- 
los which can help alleviate the constraints from gamma 
rays produced in our own galaxy. One is the possible 
presence of many subhalos with low velocity dispersion 
being the principal regions of dark matter annihilation 
would displace the gamma rays away from the galactic 
center, where constraints from HESS are strongest. The 
other arises from new studies of the effects of baryons 
on the DM velocity dispersion profiles, which imply that 
Sommerfeld-enhanced annihilation would be suppressed 
near the galactic center. These two effects could work to- 
gether such that PAMEL A/Fermi observations are dom- 
inated by subhalo annihilations, while the the 511 keV 
excess is enhanced by the larger DM velocities in the 
galactic center. We have ruled out a third possibility 
in the present class of models, namely that the interme- 
diate particles decaying into e + e~ travel away from the 
galactic center before decaying; we showed that such par- 
ticles would necessarily spoil BBN because of their long 
lifetimes and high abundance in the early universe. 

All of these loopholes are relatively unimportant for a 
related class of constraints, which considers the gamma 
rays emitted by all halos at all redshifts, including their 
effect on the CMB. It would be important to reconsider 
these constraints specifically for annihilations in which 
XX ~> 4e, and higher numbers of electrons/positrons, due 
to cascading of the dark gauge bosons, to see how much 
they really constrain the present class of models. 

We have shown that there is a new potential signal 
which could provide additional evidence for nonabelian 
DM, if the XDM interpretation of the INTEGRAL 511 
keV excess is correct. The DM transition magnetic mo- 
ment interaction induced at one loop, due to the non- 
abelian terms in the gauge kinetic mixing, should give 
rise to a narrow gamma ray line with energy equal to the 
mass splitting of the two DM states, expected to be of 
order a few MeV. Our estimates show that for a g ~ 0.04, 
the strength of this line can be close to the current sensi- 
tivity of INTEGRAL/SPI, if the phase space for excited 



dark matter decay is accidentally small. 

Another potential signal is through the couplings of 
the portal particle which interacts with the SM, with 
coupling reduced by the factor e if it is a dark gauge bo- 
son, or mixing angle 9 if it is a Higgs boson. The e factor 
is already strongly constrained, depending on the mass 
fi of the portal boson, by precision QED tests and beam 
dump experiments. Proposed fixed target experiments 
could further probe the allowed range of e in the near 
future. We have derived the bound t^4x 10~ n in the 
case of gauge kinetic mixing, from the requirement that 
dark gauge bosons annihilate to e + e~ efficiently before 
nucleosynthesis. If any of the gauge bosons are stable, we 
get the stronger bound e >, 3 x 10 -10 from overclosure. 
This is three orders of magnitude smaller than current 
limits from the E137 beam dump experiment. In the 
case of Higgs mixing, we find an analogous lower bound 
9 >, 10~ 6 on the Higgs mixing angle. 

We have also shown that Higgs mixing rather than 
gauge kinetic mixing can be a viable portal to the SM, 
despite first appearances that it would tend to induce too 
high a rate of XiXi ~> //) where / is any SM fermion, 
preferably the top quark. Such diagonal couplings of the 
Higgs to XiXi might also be expected to lead to too high 
of a cross section for scattering from nucleons, \N — > 
xN', in direct DM detectors. We have shown that in 
fact the constraints can easily be satisfied for reasonable 
values of the couplings. 

Concerning specific models, the simplest possibility is 
DM in the doublet representation, which has the poten- 
tial for implementing the XDM mechanism to explain IN- 
TEGRAL/SPI observations of 511 keV gamma rays from 
the galactic center. We noted that the doublet model 
must get its mass splitting from a Yukawa coupling to 
a triplet Higgs, rather than from radiative corrections, 
which divorces the scale of the mass splitting from the 
strength of the gauge coupling. 

A recurring theme of our paper was the possibility of an 
inverted mass hierarchy in models with three or more DM 
components, which can help boost the production of the 
511 keV excess by the XDM mechanism. A Z2 discrete 
symmetry is required to keep the intermediate-mass DM 
state stable in this scenario. We showed that models of 
triplet and quintuplet DM afford many examples with the 
desired properties, as well as the alternative normal mass 
hierarchy. In both triplet and quintuplet DM models, 
it is possible to choose Higgs VEVs such that the mass 
splittings are smaller than the generic estimate a g fi for 
the size of the radiative corrections, which can be help 
explain the hierarchy between splittings needed for the 
iDM or inverse mass hierarchy XDM scenarios. An in- 
depth reanalysis of the viability of the XDM mechanism 
for explaining the 511 keV observation, in the context of 
the present class of models, is in progress [o9| . 

We have not focused on the DAMA/LIBRA annual 
modulation as one of the signals to be explained by non- 
abelian dark matter, even though that was cited as one 
of the original motivations for this class of models. Our 
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choice stems from recent results [52| which note that 
the iDM explanation for the signal is ruled out at the 
99% c.l. by data from the ZEPLIN-II and at 95% c.l. 
by XENON10 and CRESST II observations, for DM 
whose mass is in the range of interest for explaining the 
PAMELA/Fcrmi lcpton excesses. Another reason is that 
our proposal of the inverted mass hierarchy for boosting 
the XDM mechanism is at odds with the normal mass 
hierarchy needed for iDM. However if one believes there 
is still room for iDM to work, then the desired mass split- 
tings can be achieved within the models considered here 
(for DM in triplet and higher representations), since there 
is great freedom to adjust the DM spectra through ratios 
of the VEV's of the dark Higgs fields. 

In summary, SU(2) gauge theories of DM continue to 
offer an elegant explanation for numerous effects, with 
intricate implications for cosmology, DM halo properties, 
laboratory tests, and the prospect for being ruled in or 
out in the near future. 
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APPENDIX A: TRANSITION MAGNETIC 
MOMENT 

By routing the external momenta through the loop in 
the appropriate way, the expression for the loop diagram 
which gives rise to the transition magnetic moment can 
be written as 



d 4 p 



(2tt)4 [(p + qf - M?][(p + 5qf - ^}[(p - 5qf - ^} 

(Al) 

antisymmetrized over /i, v (since it is contracted with 
F^ u of the external photon), where q is the average 4- 
momentum of the two external DM states, q — 5(93+92), 
and 5q is half the 4-momentum of the photon, Sq = 
|(?3 — 92)7 and Mi is the mass of the virtual DM parti- 
cle. We have ignored mass differences between the two 
gauge bosons in the loop since this has a subleading ef- 
fect on the result. Using Feynman parameters and Wick 
rotating, the p integral can be done, leading to 
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dx I dy 
Jo 

in({% + y)j +(y- x )&j + Mi) 7^ 
z 2 M 2 + \{y- x)z5M$ 2 + 5Mf 23 z + ^ 2 (x + y) 



(A2) 



y, 6M* 23 = Ml - |(M| + Ml), 



where z — 1 — x — y, SM'{ 23 = M\ 
and <5Mf 2 = Mf — Mf. For the parameter values of 
interest, we find that it is a good approximation to set 
&M32 = SM 2 23 = in the denominator. By anticom- 
muting gamma matrices in the numerator and using the 
Dirac equation for the external spinors, one can show 
that § -> -|(M 2 + M 3 ), while 5f gives a subleading in 
SM23 contribution which can be neglected. Furthermore, 
it is a good approximation to set (x + y) — 1 for the co- 
efficient of /x 2 in the denominator. In this way one can 
get the analytic approximation which we have nu- 
merically verified to be good in the range of parameters 
of interest. 



APPENDIX B: RADIATIVE DM MASS 
SPLITTINGS 

In this appendix we present the radiative mass cor- 
rections to a DM multiplet Xi by virtual massive gauge 
bosons, as shown in fig. 03(a). Assume the DM multiplet 
transforms under a gauge group G with generators Tj. 
Also assume the gauge bosons Aj have mass fij. Figure 
02a) gives a correction to the self-energy of \% °f 



SMi = g 



2 Jij rpj 

/ j ia ai 



dx 



• tE 4M a - 2p 



(27T) 



A) 2 



(Bl) 



where A = -M 2 (l-a;)a;+^ 2 a;+M 2 (l-x) ~ M 2 (l-a;) 2 + 
n'jX. After integrating over the Euclidean 4-momentum 
k,E and using the equation of motion to set i> — > Mi, we 
find two pieces, one of which is ultraviolet (UV) diver- 
gent, while the other is infrared (IR) divergent as ju,j — > 0, 



SMi 



9_ 

8tt 



2 E T - T « / dx ( 2M - - M i x ) ( ln A2 - ln A ) 

j,a Jo 



(B2) 



where A is the ultraviolet cutoff of kg. We are only inter- 
ested in the IR divergent term because the UV divergent 
term cancels out when considering mass splittings. No- 
tice that the IR divergence occurs when x — > 1, so we 
can set x = 1 in the first factor of the integrand. 

Further assuming that the gauge boson mass is much 
smaller than the DM mass, fij <C M a , and ignoring the \ 
mass splittings on the r.h.s., the IR divergent term turns 
out to be 

el 



2tt 



Xtf 



(1 - a;) 2 M 2 



rp 'J rp 'J 

ia ai 



+ 21n((l-x)M x ) 



(B3) 



We dropped the last term because it is the same for each 
component of the DM multiplet and thus has no effect 
on the mass splitting. 
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Now we turn to the mass correction from virtual Higgs 
bosons. For illustration, suppose we have triplet DM Xi 
coupled to a quintuplet scalar Sij, via h s XiSijXj- This 
induces radiative mass corrections through the diagram 
of fig. [5]Jb). Similarly to the gauge boson case we have, 



matrix element is 



SMi = 




+ 2\a((l-x)M j )J 
(B4) 



where the last approximation applies if the dark Higgs 
boson mass is much smaller than the dark matter mass, 
my < M a . 



APPENDIX C: 



DM ANNIHILATION CROSS 
SECTION 



Here we derive the annihilation cross sections for DM 
in the doublet, triplet and quintuplet representations. 
Because of its relative simplicity we start with the triplet 
case, assuming Majorana DM, whose gauge interaction 
Lagrangian is 



1 



' g tabc XaB^X 



(CI) 



For the annihilation channel X1X1 
interaction is 

3Xi^2 7^X3 



B2B2, the relevant 



(C2) 



where the antisymmetric property of the Majorana vector 
current, x^l^Xi — — Xi7 M X3 nas been used. The matrix 
element for X1X1 ~ y B2B2 is then 



M11 



-ig 2 v{pi)t{qi 



1 



~ig 2 v(pi)^(q 2 )- —fiq 1 )u{p 2 ) (03) 

where pi are the incoming x momenta and qi are the 
outgoing B momenta, and the e's are polarization vec- 
tors. This has algebraically the same form as the ma- 
trix element for electron-positron annihilation. In the 
nonrclativistic and M x 3> /i limits, the spin-averaged 
squared matrix element is 2g A for both the t and u chan- 
nels squared. The interference term vanishes in these 
limits. The sum over final states includes X1X1 — > B3B3 
as well, so the total is 8g 4 . 

Next we consider the X1X2 — > B\B 2 , which procedes by 
a sum of i-channel mediated by internal X3 an d s-channel 
mediated by B3 (using the 3-gauge boson vertex). The 



Mi 



-ig 2 v(p 1 ){(q 1 )- 



f(q 2 )u(p 2 ) 



. 2 y(pihuu(p 2 ) r . . 
- ig — 2 2 — ivvxips + q 2 )^ 

-ilv^ijPs+qijx + Tq^xiqi-q^u] (C4) 

where p s — Pi + p 2 - We find that the spin- 
averaged squared matrix element gets contributions of 
2<? 4 , — (19/4)g 4 and 4g 4 from the t 2 , s 2 and interference 
channels, respectively, in the same limits as mentioned 
above. (The fact that the direct s 2 term is negative is 
due to the unphysical polarizations in the sum over final 
state gauge bosons; only the full amplitude squared is 
physically meaningful.) The total is thus (5/4)<? 4 . 

Finally we must average over the initial colors to give 
(|M 2 |) = \Ml, + \M\ 2 = [8/3 + (2/3)(5/4)] ff 4 = 
(7/2)g 4 . This must be multiplied by an additional factor 
of 1/2 for the indistinguishability of the final states. The 
differential cross section is thus 



da 



(7/4)ff 4 
64nsM 2 v 2 



(C5) 



evaluated in the center of mass frame. Integrating over 
the range of the Mandelstam t variable St = AM 2 v and 
using the relative velocity v re i = 2i>, we obtain the cross 
section (|34"|) for the triplet case. 

To find the annihilation cross section for DM in other 
representations, we work out the group theory factors 
for the general case. For the t 2 and u 2 terms, these 
take the form (l/d%) X; Q6 tr(T a T a T b T fc ) = C 2 (R) 2 /d R , 
where d R = 2j + 1 is the dimension of the spin j 
representation and C 2 (R) = j(J + 1) is the quadratic 
Casimir invariant. The factors (1/dfj) come from av- 
eraging over the colors of the initial states. For the 
s 2 term we get (1/4) £y £ a6cd ^ahd^dy = 

2/d 2 R J2 c tr[T c T c ] = 2C 2 {R)/d R . For the st term, we 
find (l/4)tr(T Q T b T c )e abc = (i/2d 2 B )e abd e ab MT d T c ) ~ 
C 2 (R)/d R . 

Using these results, we can generalize the triplet com- 



putation to other representations. Using t 



and 



st to represent the contribution to the squared matrix el- 
ement from any definite external states (hence encoding 
the spin sums but not the group theory factors) we have 

(\M 2 \) = l(^(t> + u>) + S 2 + s?) (06) 

in the triplet representation, which based on the above 
group theory factors, generalizes to 



(07) 




where t 2 = u 2 — 2g 4 , s 2 
obtain 



-(19/4) ff 4 and st = 4 5 4 . We 



j(j + 1) - J 



(08) 
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for the spin-j representation. This must still be multi- 
plied by the symmetry factor 1/2 for identical final states. 



where sg = sin 6*, eg — cos 6, and 1 is the 2x2 unit ma- 
trix Treating e as a perturbation, one can perform a 3x3 
rotation to get rid of the off-diagonal e blocks, using 



APPENDIX D: RATE OF 3B -> 2B 
DOWN-SCATTERING PROCESS 



f 1 - -qrf/2 -T] 
1 77* 1-77*77/2 



(E4) 



In the case of e -C 1 where dark gauge bosons do not 
stay in kinetic equilibrium with the SM particles, the pro- 
cesses which can delay them from dominating the energy 
density of the universe after becoming nonrelativistic are 
those which convert 3 to 2 particles. The squared matrix 
element for ZB — > 2B in nonabelian gauge theory is of 
order \Ai\ 2 ~ g 6 /p 2 if the B's are nonrelativistic. The 
Boltzmann equation for the number density of B's takes 
the form fig + 3HriB = C. The collision term of interest 
for 3B -> 2B is 



C 



n 



(2tt) 3 2£, 

(2tt)^( Pi + P2 + P3 
n|(r)(3/2) 5 (47r) 4 g 6 
(27t) 11 /j 3 P 2 
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where the nonrelativistic density is ns{T) = 
(/xT) 3 / 2 e -AI / T . To convert this to a rate, we should 
divide by one power of tlb- Equating this to the Hubble 
rate, we find that 



2x f 



In 



3.4a 3 M p 



-In Xf 



(D2) 



leading to the result in the text above eq. (|40p . 



where 77 = (A — B x l) _1 e = A _1 e. Under this rotation, 
the 2x2 block A transforms again, receiving a correction 
A — > A + 6 A of the form 



5 A = ~\{A, A^XA- 1 } + {X, A- 1 } 



BAXA~ 



_(E5) 

where X is the 2x2 matrix ee* , whose only nonvanishing 
component is e 2 in the 2,2 position. We find that to 
leading order in S, 



5 A = 







c e sgSe- 



A-B 



2(A-B) 2 



1 - 



2(A-B) 



(E6) 



where A — g v . B gets a similar correction, but the 
correction to A is more important since this splits the 
gauge boson mass eigenvalues, whereas B is already well 
separated from the eigenvalues of the A matrix. The 
final step for the gauge boson mass eigenvalues is to di- 
agonalize A + 5 A. The off-diagonal elements of 8 A give 
only 0(S 2 ) corrections to the gauge boson mass split- 
tings, which we are ignoring; thus, denoting A = g 2 v 2 



and B — y v 



y ■ 



the resulting masses are 



APPENDIX E: DIAGONALIZATION OF 
SU(2)xU(l) MODEL MASS MATRICES 

We give the details of approximately solving for the 
mass eigenvalues and eigenstates in the SU(2) x U(l) DM 
model discussed in section llXl The dark gauge bosons are 
denoted by Bi and Y, and the effects of mixing with the 
SM hypercharge are neglected in the following approx- 
imations. Since Bi does not mix with the other fields, 
we can consider the nontrivial 3x3 mass matrix in the 
B\, B3, Y basis, writing it in the form 



where A = di&g(g 2 v 2 , g 2 v 2 



A e 
e 4 B 
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S) 



gy 
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'■i> \ ' • 

y"v" + y l2 4> 2 ■ It is convenient to change bases 
using a global SU(2) rotation which makes only the lower 
component of e nonzero. This is accomplished using a 
2x2 rotation R(6) in the B1-B3 subspace, with tan# = 
t\jt2- In the new basis, A is no longer diagonal, 

\ - ,,-r-i + ,) [ s « ! (E3) 
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f4 
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A + s 2 6 + 0(S 2 ) 
A + S 

A+^-^ + 0(6) 
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(E7) 
(E8) 

(E9) 

(E10) 



The relation between the flavor (unprimed) and mass 
(primed) eigenstates is B2 = B' 2 and 



(Ell) 









-( 








where 



A + 5-B : 



i< = CgSg(A — B) 



(E12) 



Here ip is the small angle of the rotation which diagonal- 
izes (|E6|) . It gives rise to 0(5 2 /e 4 ) contributions to the 
mass splittings of the DM states, which are larger than 
the 0(S 2 ) terms we have ignored thus far. These formu- 
las assume 6 < e 2 /\A — B\, so they are not valid in the 
limit e — > 0. 
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